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Abstract 

This paper focuses on the impact of Racah on crystal- and ligand-field theories, 
two branches of molecular physics and condensed matter physics (dealing with ions 
embedded in aggregates of finite symmetry). The role of Racah and some of his stu- 
dents in developing a symmetry-adapted weak-field model for crystal-field theory is 
examined. Then, we discuss the extension of this model to a generalized symmetry- 
adapted weak-field model for ligand-field theory. Symmetry considerations via the 
use of the Wigner- Racah algebra for chains of type SU (2) D G is essential for these 
weak-field models. Therefore, the basic ingredients for the Wigner-Racah algebra 
of a finite or compact group are reviewed with a special attention paid to the SU (2) 
group in a SU{2) D G basis. Finally, as an unexpected application of nonstandard 
SU(2) bases, it is shown how SU(2) bases adapted to the cyclic group allow to build 
bases of relevance in quantum information. 



1 Introduction 

The legacy of Giulio Racah (Firenze, 1909-1965) stems mainly from his four papers Theory 
of complex spectra published between 1942 and 1949 [H [21 El H], his notes on group- 
theoretical methods in spectroscopy based on lectures given at the Institute for Advanced 
Study in Princeton in 1951 [6], and his book on irreducible tensorial sets written in 
collaboration with his cousin Ugo Fano [TJ. 

It was the purpose of the first paper of his famous series [U 12 13, 4J to substitute to the 
numerical methods of Slater, Condon and Shortley general methods more conformable to 



the Dirac representation of state vectors. The main achievements realized in [H [2J |3l HJ 
El [7] deal with irreducible tensor methods, Wigner-Racah algebra (or Racah-Wigner 
algebra or Racah algebra, a concept to be precisely defined in Section 3) and group- 
theoretical methods involving chains of Lie groups. More precisely, let us mention the 
following important contributions. 

• The development of the algebra of coupling and recoupling coefficients for the SU(2) 
group in a SU{2) D U(l) basis, with introduction of the V and V functions (the V 
symbol is identical to the 3—jm Wigner symbol up to a permutation of its columns) 
and of the W, W and X functions (the W and X symbols are identical to the 6-j 
and 9-j Wigner symbols, respectively). 

• The introduction of the concept of a SU(2) irreducible tensor operator that gener- 
alizes the notion of a vector operator and the generalization to tensor operators of 
the Wigner- Eckart theorem for vector operators. 

• The introduction of the notion of a unit tensor operator, the matrix elements of 
which in a SU(2) D U(l) basis are nothing but Clebsch-Gordan coefficients (up to 
a multiplicative factor), with the advantage that any tensor operator is proportional 
to a unit tensor operator. 

• The introduction of the concept of seniority which is related to the state labeling 
problem. 

• The development of the notion of coefficients of fractional parentage, previously 
introduced by Goudsmit and Bacher, which make it possible to develop a n-particle 
wavefunction in terms of (n — l)-particle wavefunctions. 

• The introduction of chains of Lie groups, involving both invariance and classification 
groups, for characterizing state vectors and interactions involved in spectroscopic 
problems. To implement the use of chains of groups, Racah introduced a factor- 
ization lemma and developed the notion of a complete set of commuting operators 
(involving Cartan operators, invariant or Casimir operators and labeling operators) 
in a group-theoretical context. 

The series of seminal works [U [21 El El El El E] opened the way for many applications 
by Racah himself, his students and a large part of the community of scientists working in 
atomic and nuclear spectroscopy (see the list of Racah's publications in [8]). In particular, 
the methods of Racah were popularized by Judd [9J, Wybourne [ID] , and Condon and 
Odaba§i [11] in atomic physics and by de-Shalit and Talmi [12] in nuclear physics (see 
also |13[ dH [T51 [TB] for recent developments in nuclear and molecular physics). The 



basic concepts introduced and/or developed by Racah in his pioneer works were also of 
considerable importance in molecular and condensed matter physics. More specifically, 
these works stimulated an enormous quantity of developments in vibration-rotational 
spectroscopy of molecules and in crystal- and ligand-field theories as will be shown below. 

We shall be concerned here with the impact of Racah on crystal- and ligand-field 
theories, two theories which deal with optical, magnetic and thermal properties of ions 
embedded in molecular, condensed matter or biological surroundings. Racah never pub- 
lished papers about these theories. However, he was interested in molecular physics as 
shown by the fact that he published in 1943 a short note on the structure of the Mo(CN)g~ 
complex ion [T7]. His interest for molecular physics and the physics of ions in crystals was 
reinforced and stimulated by a seminar given by his colleague Willy Low in the Depart- 
ment of Physics of the Hebrew University of Jerusalem in 1956 [T81 [T9"} ¥2U\ . The seminar 
was devoted to the role of crystalline fields on the optical spectra of transition-metal ions 
(like Ni 2+ and Co 2+ ) in crystals. Racah became very much interested in this research 
subject and decided to guide students in this direction. His idea was to combine his ir- 
reducible tensor methods with the group-theoretical methods largely used in crystal-field 
theory (but limited in those times to a qualitative explanation of the level splitting for a 
given ion embedded in a finite symmetry surrounding). Along this line, Racah and Low 
directed two graduate students, Schoenfeld who studied the case of the d 2 and d 3 config- 
urations in cubic symmetry [2T] and Rosengarten who dealt with the case of d A and d 5 
configurations in the same symmetry [22]. Then, Racah asked another student, Flato, to 
work out the more involved case of the d 2 and d 3 configurations in trigonal and tetragonal 
symmetries [23J. Five years after having completed his thesis, the material contained in 
Flato's thesis was still of such an interest that he was asked to publish it (for the main 
part) [2l] (see also [25J ) . Research in that direction continued with a general formalism 
and a symmetry-adapted weak-field model developed by the present author in his thesis 
prepared under the guidance of Flato [26, EHj- 

It is one of the aims of the present review to show how Racah directly and indirectly 
contributed to the penetration in crystal- and ligand-field theories of the tools he originally 
developed for atomic and nuclear spectroscopy. Another aim of this article is to show 
how the Wigner-Racah algebra for a group of molecular or crystallographic interest can 
be deduced from the one of SU(2) in a nonstandard basis. To a large extent, this paper 
constitutes a brief review of the methods and models used in crystal- and ligand-field 
theories as well as a pedestrian presentation of the Wigner-Racah algebra for a chain 
of groups involving finite and/or compact groups. As an application of the SU(2) D G 
chain, where G is a cyclic group, a brief contact is established with quantum information, 
a field of considerable interest in the present days. 

The material in this paper is organized as follows. Section 2 deals with crystal- and 



ligand-field theories. The basic ingredients for the Wigner-Racah algebra of a finite or 
compact group together with some illustrative examples are given in Section 3. Section 
4 is devoted to a short incursion in quantum information via the use of specific chains of 
type SU(2) D G. 

Most of the notations are standard. The star denotes complex conjugaison, 5^ the 
Kronecker delta symbol of a and b, and A' the adjoint of the operator A. We use a 
notation of type (as in Racah's papers), or \ip) (as in quantum information), for a 
vector in an Hilbert space and we denote (<f)\if)) and \<p)(ip\ respectively the inner and 
outer products of the vectors and \<p). Specific notations on group theory shall be 
introduced later. 

2 Crystal- and ligand-field theories 

2. 1 Generalities 

Crystal- and ligand-field theories deal with the description and interpretation of electronic 
and magnetic properties (optical spectra, electron paramagnetic resonance spectra, photo- 
electron spectra, etc.) of a partly-filled shell ion in a molecular, condensed matter or 
biological environment. Crystal- field theory (based on the use of atomic orbitals) goes 
back to the end of the twenties with the seminal work by Bethe [29J and was applied to 
the spectroscopy of ions in solids in the early days of quantum mechanics. It is only in the 
fiftees that ligand-field theory (based on the use of molecular orbitals) was the object of 
numerous studies. In modern parlance, crystal- and ligand-field theories are special cases 
of the theory of level splitting. 

As a typical example, let us consider the case of a ruby crystal. It consists of corindon 
(AI2O3) doped with trivalent chromium ions (Cr 3+ ) in substitution with trivalent alu- 
minium ions (Al 3+ ). The electrons of each Cr 3+ ion are thus subjected to inhomogeneous 
electric fields arising from the ligands or coordinats constituted by the oxygen atoms. 
These electric (or crystalline) fields yield a level splitting of the energy levels of the Cr 3+ 
ion. One-photon transitions in the visible between the split levels are responsible for the 
nice pink to blood-red color of ruby. 

The distinction between crystal-field theory and ligand-field theory is as follows. In 
crystal-field theory one uses atomic orbitals for the central partly-filled shell ion (the Cr 3+ 
ion in our example) whereas in ligand-field theory one considers molecular orbitals made 
of linear combinations of atomic orbitals of the central ion and of the ligands or coordinats 
(the O 2- ions in our example). 



2.2 The Hamiltonian 



We shall consider the common case of an ion with a £ N atomic configuration (N equivalent 
electrons on a nC. shell outside of a set of closed shells). The I = 2 case corresponds 
to transition metal-ions and the I = 3 case to rare earth and actinide ions. In first 
approximation, the perturbation Hamiltonian H for such an ion embedded in a crystalline 
field reads 

H := He + Hso + H c f (1) 

where He stands for the two-body Coulomb interaction between the N electrons, H so 
the one-body spin-orbit interaction for the N electrons and H c f the one-body interaction 
between the N electrons and the environment of the central ion. Obviously, He an d 
H so are rotationally invariant and H c f is invariant under the point symmetry group G 
of the ion and its surrounding. Therefore, the 0(3) dG chain of groups naturally plays 
an important role in the description of the ion in its environment (the three-dimensional 
orthogonal group 0(3) is isomorphic with the three-dimensional rotation group). When 
G contains only rotations, it is sufficient to consider the SO (3) D G chain (SO (3) is 
the three-dimensional special orthogonal group) or the SU(2) D G* chain, where SU(2) 
and G* are the spinor groups (double groups in the terminology of Bethe) of 50(3) ~ 
SU(2)/Z 2 and G ~ G*/Z 2 , respectively. We can thus understand the importance of both 
continuous and finite groups in crystal- and ligand-field theories. 

In view of the various terms in H, we can have several families of models. The 
situations 

He > H so > H c f (2) 

and 

H c f > He > H so (3) 

correspond to the so-called weak-field model and the strong-field model, respectively. The 
strong-field model was mainly developed in the fifties by Tanabe, Sugano and Kamimura 
in Japan [30j SB [321 [331 HQ ES] and by Griffith in England [361 EH EHJ EH EH SB S21 El 
H5] , and later by Tang Au-chin and his collaborators in China El SSI SU ED] as well as 
by Smirnov and his collaborators in the former USSR [SB E21 ESI EU ES]- The weak-field 
model, although worked out in the early days of crystal-field theory, was systematically 
developed from the sixties. In particular, a symmetry-adapted version of the weak-field 
model was introduced, as we mentioned in the introduction, following a suggestion of 
Racah by two of his students, Schoenfeld [21] and Flato [2S1 [21] • It was further developed 
by the present author and some of his collaborators [261 I2H Ell EB ESI E9] (see also 



HS1 STl SH S9l [50] and [5U E21 [531 EH E5] ) . In crystal-field theory, the weak- and strong- 
field models are a priori equivalent if the matrix of H is set up on the 
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state vectors of the £ N configuration. Nevertheless, the implementations of the two models 
are quite different as it will be shown below. As an illustration, we shall now discuss in 
turn the two models (strong- and weak-field models) in the special case of d ions in cubic 
symmetry. 

2.3 Strong-field models 

It is difficult to describe the strong- field model in the general case of £ N in G. Hence, we 
consider the case of a d ion (£ = 2) in octahedral symmetry (G = O). The restriction 
50(3) — > O yields the following decomposition 

2 = E®T 2 (5) 

of the irreducible representation class (IRC) of 5*0(3) associated with 1 = 2 into a direct 
sum of the IRCs E and T 2 of finite group O. As a consequence, there is a splitting level: 
the five degenerate d orbitals give rise to a E doublet (with two degenerate orbitals e) and 
a T 2 triplet (with three degenerate orbitals t 2 ). The e and t 2 orbitals can be considered 
as symmetry-adapted atomic orbitals (in crystal-field theory) or as molecular orbitals (in 
ligand-field theory). The distribution of the N electrons on the t 2 and e orbitals, according 
to the Pauli exclusion principle, yields (molecular) configurations of type t^e x . Then, 
we can form (molecular) terms t2(SiTi)e N ~ x (S 2 T 2 ), where Si and 5*2 are the total spins 
for the x and N — x electrons on the t 2 and e orbitals, respectively. Furthermore, Ti 
(contained in T® x ) and Y 2 (contained in E®^ N ~ xS> ) denote the IRCs characterizing the 
orbital parts of the t 2 and e electrons. The next step is to couple 5*1 with 5*2 to get the 
total spin S (contained in S\ <S> S 2 ) and Ti with T 2 to obtain T (contained in Ti ®r 2 ). This 
leads to (molecular) states t|(S'iri)e Af ~ a: (S'2r2)ST. Finally, the coupling of S (decomposed 
into IRCs of O*) with T gives the total IRC Tt (an internal branching multiplicity label 
b is necessary when occurs several times in the reduction of S <E> T). As a result, we 
get state vectors of type 

|^(5 1 r 1 )e Af ^(5 2 r 2 )sr6r T7T ) (6) 



which are expressed (via complicated formulas) in terms of one-electron state vectors by 
means of coupling coefficients and coefficients of fractional parentage. Note that the label 
7t i n © is necessary when the dimension of is greater than 1. 



The calculation of the matrix elements of % c f in the strong-field basis ([6]) is elementary. 
However, this is not the case for % c + H so - The construction of the matrix of He + H S o 
on C^ e+2 state vectors (jHJ) requires the knowledge of coupling and recoupling coefficients 
for both 577(2) and G* as well as coefficients of fractional parentage for the configurations 

From the practical point of view, the just described strong-field approach leads to: 

• a five-parameter model in a crystal-field framework where the t% and e orbitals are 
atomic orbitals, called ordinary strong- field model, with 3 parameters for T-Lq (F , 
F 2 and F 4 of Slater or A, B and C of Racah, see the appendix), 1 parameter for 
Hso (Cnd) an d 1 parameter for % c f (10/Jg) 

• a fourteen-parameter model in a ligand-field framework where the t% and e orbitals 
are molecular orbitals, called generalized strong-field model, with 10 parameters for 
He, 2 parameters for % so and 2 parameters for % c f. 

The strong-field models present several drawbacks. The case of d in O is difficult to 
extend to the case of £ in G: replacing O by G and/or d N by £ requires that the calcu- 
lation for Tic and H so , which involves complicated Wigner-Racah algebra developments 
for the G or G* group with several phase problems, have to be done again. This kind of 
difficulty does not appear in a weak-field approach as shown below. 

2-4 Weak- field models 

In the case of £ in G, we may think to use atomic state vectors of type \n£ N aSLJ M). 
However, such state vectors, adapted to the SU(2) D U(l) chain, are not generally 
adapted to the G* symmetry group. The idea of Racah was to use linear combinations 
of the vectors \n£ N aSLJ M) transforming as IRCs of G* and to employ his methods 
for calculating the energy matrix of H. Therefore, the matrices for Tic and H S o, in a 
577(2) D G* symmetry-adapted basis, are the same as the ones of atomic spectroscopy 
(already calculated by Racah or easily calculable from Racah's methods) and the matrix of 
H c f depends on reduced matrix elements of one-electron Racah unit tensor operators and 
SU(2) D G* symmetry-adapted Clebsch-Gordan coefficients. Thus, the implementation of 
the symmetry-adapted weak-field model is easier than the one of the ordinary strong-field 
model. Following Racah's idea, Schoenfeld and Flato calculated the matrix of H for the 
d 2 and d 3 configurations in cubic symmetry [21] and in tetragonal and trigonal symmetries 
[23| |23] . Later, Low and Rosengarten dealt with the case of the d 5 configuration in cubic 
symmetry in connection with the optical spectra of Mn^ and Fe 3+ ions in crystalline 
fields [22]. 



The Wigner-Racah algebra for the 577(2) group in a SU (2) D G* symmetry-adapted 
basis of interest for the symmetry-adapted weak-field model was developed by the present 
author (26] E7J EE] [60] and further considered by several authors [EEJ E2J [63] E3 ESI ES 
1671 168] [69] . The main ingredients of the resulting symmetry-adapted weak-field model for 
£ in G can be summed up as follows. 

The symmetry-adapted weak-field state vectors are of type 

j 

\n£ N aSLJaT-i) := \n£ N aSLJM)(JM\JaT 1 ) (7) 

M=-J 

where V is an IRC of G*, a a branching multiplicity label to be used when the (J) 
IRC of SU(2), associated with the J quantum number, contains T several times and 
7 a multiplicity label to be used when the dimension of the T IRC is greater than 1. 
In (J7J), the (J M\JaT^) reduction coefficients are elements of a unitary matrix which 
reduces the representation matrix associated with the (J) IRC of SU{2) into a direct 
sum of representation matrices of G*. They have to be distinguished from the reduction 
coefficients obtained from the diagonalization of an operator invariant under the G group 
[70] ITT] [721 175] ITT] [75] . The (JM\ JaTj) reduction coefficients are chosen in such a way 
that the set 

{\n£ N aSLJaT"f) : 7 ranging} (8) 

spans a representation matrix associated with T independent of the atomic quantum 
numbers and that the values of the corresponding coupling coefficients (the / coefficients 
below) are square roots of rational numbers. Then, the matrices for He an d "H S o follow 
from 

(n^aSLJaY^Uc^a'S'L'J'a'Y 1 ^) = 5ss'&Lvhj>&< M >&m& ri ' 

xA(S, L, J){nl N aSLM s M L \U c \nl N a' SLM S M L ) (9) 

and 

(n£ N 'aSLJaTj\n ao \ni N 'a' 'S' ' V 'J'a'T'i) = 5 J , 7 ><W<W T / 

x(n£ N aSLJM\H so \n£ N a'S'L'JM) (10) 

where A(S, L, J) is 1 if S, L and J satisfy the triangular condition and otherwise ; in 
OH]) and f llOp . the matrix elements in the right-hand sides are independent of the magnetic 
quantum numbers Ms, Ml and M, respectively. Clearly, the energy matrices for Tic and 
l-L so do not depend on the G group and are easily builded from the works of Racah (the 
matrix elements in the right-hand sides of OH]) and ([TO]) are known for the p N , d N and / 



configurations [76] or easily calculable from computer programs). On the other hand, the 
matrix of T-i c f can be readily set up by making use of the development 

n cf = J2D[ka ]U^ o (11) 

kao 

where is a component of a Racah unit tensor operator XJ k invariant under G (i.e., 

transforming as the r identity IRC of G). In (TTT]) . D[ka ] are crystal-field parameters 
connected to the B k parameters (in Wybourne's normalization (HUE?]) via 

D[ka ] = (-l) e (2e+l) J ^ B k q (kq\ka T olo y (12) 

q= — k 

and ao is a branching multiplicity label to be used when T appears several times in the 
decomposition of the (k) IRC of SO(3). (The index 7 in ffTTj) and (fT2l) is not really 
necessary since T is a one-dimensional IRC; it is mentioned only for aesthetic reasons.) 
Then, the matrix elements of "H c / in a SU(2) D G* symmetry-adapted weak-field basis 
are given by 

(n^aSLJaT^Uc^n^a'S'L'J'a'T'i) = <W<W 7 y 
x (-1) S+L '+ V(2J + 1)(2 J' + 1) ^ D[A;a ] (13) 

fcao 

x (^aSLH^H^a'SL') ( ^ ^ ^ 1 / T ^ J * 
V " " \j' S J j \aT a'T a T 

where {• • •} stands for a 6-j Wigner symbol and / is a coupling coefficient defined by 

f( J J k ) := Yl E X)(-/M|Jor 7 )*(*9|*aor 7o)(J / W|J / a , r7) 
\ai at a l / m=-j m'=-j> g =-k 

T / J k J' \ 

x (-1) J (14) 
V-M g MV 

This / coefficient is independent of 7 [211 [26] . It is a particular case of the / coefficient 
defined in [26J by 

(j j k\ ^ ^ 2 k 

J := 52 52 {jim 1 \jmT{kq\kfJ 1 ){j 2 m2\j 2 fJ'2) 
r r r ' mi=— ji m.2=— J2 q=—k 

x ( Jl fc j2 ^ (15) 

V -mi q m 2 J 

where 

Mi := 01^7!, /i 2 := a 2 r 2 72, M := ar 7 (16) 



(see also [60J). As a conclusion, the calculation of the matrix of T-L in a symmetry- adapted 
weak-field basis (via <Q, ( ITU]) and ( ITBl ) is considerably simpler than in a strong- field basis. 



In the case of d in O, the just described symmetry- adapted weak- field approach, 
based on (Q, ffTO]) and (fT3l) . leads to a weak- field model which is equivalent to the ordinary 
(or ionic) strong-field model with the parameters A, B, C, ( n d and lODq. More generally 
for £ in G, the symmetry-adapted weak-field model and the ionic strong-field model 
are equivalent. However for £ N in G, the symmetry- adapted weak- field model and the 
generalized (or covalent) strong-field model are not equivalent. Thus, it is desirable to 
develop a generalized symmetry-adapted weak-field model equivalent to the generalized 
strong-field model. This will be done in the next section. 

2.5 Generalized weak-field model 

To generalize the symmetry-adapted weak-field model, we keep the symmetry-adapted 
weak-field basis ([7]) intact in order to take advantage of its simplicity. The sole modifi- 
cation to be done consists in replacing the Hamiltonian H, by an effective Hamiltonian 
H e ff. The Hamiltonian H c r for £ in G should reduce to "H for some special values of its 
parameters, should be an Hermitian operator invariant under the G group and the time- 
reversal operator, and should contain one- and two-body spin and orbit interactions. For 
the sake of easy calculations, H e g should involve a coupling scheme which is reminiscent 
of the {SLJ} coupling scheme of the state vectors ([7]). Therefore, we take H e n in the 
form 

H cS := D[(kik2)k s {k3h)k L ka ) 

i,j all k ao 

x {{u™(i) ® u^(j)} te) <8> {u<*»>(») ® u^) (j )}(M } « o7o (17) 

where the u's are one-electron Racah unit tensor operators with {u^ kl \i) £g> u^ 2 ^ (j)}*- fcs '- ) 
acting on the spin part and {u^ 3 ^(i) ®u( fc4 )(j)}( fei ) on the orbital part of the state vectors 
(JTJ) . The sums over i and j in ffTTj) are extended over the N electrons and the sums 
over the fc's and ao are limited, like in (TIT]) , by hermiticity and symmetry properties of 
H c g (invariance under the G group and the time-reversal operator) and by the selec- 
tion rules on the matrix elements of H e s in the basis ([7]). Furthermore, the parameters 
D[(k 1 k2)ks(k 3 k4)kLkao] comprise the Coulomb interelectronic parameters, the spin-orbit 
parameters and the crystal-field parameters of the ordinary weak-field model plus some 
additional parameters to be described below. The most important (as far as a comparison 
with the generalized strong-field model is in order) parameters in H e ^ can be classified in 
the following way. 

1. The D[(00)0(/cfc)00] parameters correspond to the ordinary or isotropic Coulomb 
interaction between the N electrons. 



2. The D[{QQ)Q{k^,ki)kLkLao\ parameters with ki ^ correspond to anisotropic Coulomb 
interactions between the N electrons or ligand-field correlated Coulomb interactions. 

3. The D[(ss)l(££)10] parameter corresponds to the ordinary or isotropic spin-orbit 
interaction for the N electrons (s = 1/2). 

4. The D[(ss)l(tt)ki,kaQ\ parameters with fc^ ^ 1 correspond to anisotropic spin-orbit 
interactions for the N electrons or ligand-field correlated spin-orbit interactions. 

5. The D[(ss)0(££)kLkL(io\ parameters correspond to the ligand-field interaction. 

The building of the energy matrix of H e s in the basis (JTj) is very simple. Indeed, we 
have the following matrix elements 

(n£ N aSLJaTj\H eS \ne N ' ol S' U J'a'Y'i) = 5 T t'S 77 > 

x ^^C[(W 2 )fc s WWao]/( J J k r )H ( 18 ) 

x (n^aSLJUiu^ii) ® u (fc2) (j)} (%) ® {u^(i) g) u {ki \j)} {kL) } {k) \\n£ N a' S' L 1 J 1 ) 

where the reduced matrix element (|| ■ ■ ■ ||) can be calculated from the Racah's standard 
methods. 

The symmetry-adapted weak-field approach based on f|T7|) and (TTH|) leads to a model 
that turns out to be equivalent to the generalized strong-field model. However, the gener- 
alized symmetry-adapted weak-field model contains more parameters than the generalized 
strong- field model (e.g., the Hamiltonian given by (1T7|) contains spin-spin and orbit-orbit 
interaction parameters that do not occur in the generalized strong-field model). The 
D[(kik2)ks(k 3 k4)kLkao\ parameters can be considered as phenomenological global param- 
eters to be fitted on experimental data. All or part of these parameters can be interpreted 
and calculated in the framework of ab initio microscopic models as for instance the an- 
gular overlap model [TH1 ES], the superposition model [SD] and the MO-LCAO model 
[8T| |82| [83j El]. (See the appendix for the connection between the isotropic Coulomb 
interaction parameters and the Slater-Condon-Shortley parameters.) Of course, the gen- 
eralized symmetry-adapted weak-field model gives back the ordinary symmetry-adapted 
weak-field model as a particular case when some parameters vanish. 

By way of illustration, let us consider the case of d N in O. The corresponding Hamil- 
tonian H e g can be restricted to an operator containing 14 parameters, namely, 

• 10 Coulomb parameters: 

D[(00)0(00)00], ^[(00)0(22)00], D [(00) 0(44) 00], 
^[(00)0(04)44], ^[(00)0(22)44], £[(00)0(24)44], 
£[(00)0(44)44], £[(00)0(24)66], £[(00)0(44)66], £[(00)0(44)88]; 



2 spin-orbit parameters: 
£[(ss)l(22)10], £[(ss)l(22)34]; 



• 2 ligand-field parameters: 
£[(ss)0(22)00], £[(ss)0(22)44]; 

It can be shown that the generalized symmetry-adapted weak-field model with these 14 
parameters is equivalent to the generalized strong-field model for d N in O [58]. Such 
an equivalence was also worked out for the case of / in O. In this case, the gener- 
alized symmetry-adapted weak-field model can be restricted to involve the following 33 
parameters 

• 26 Coulomb parameters: 

£ [(00) 0(00) 00], £[(00)0(22)00], £[(00)0(44)00], £[(00)0(66)00], 
£[(00)0(04)44], £[(00)0(22)44], £[(00)0(24)44], £[(00)0(26)44], 
£[(00)0(44)44], £[(00)0(46)44], £[(00)0(66)44], £[(00)0(06)66], 
£[(00)0(24)66], £[(00)0(26)66], £[(00)0(44)66], £[(00)0(46)66], 
£[(00)0(66)66], £[(00)0(26)88], £[(00)0(44)88], £[(00)0(46)88], 
£[(00)0(66)88], £[(00)0(46)99], £[(00)0(46)10,10], £[(00)0(66)10,10], 
£[(00)0(66)12, 12a], £[(00)0(66)12, 126]; 

• 4 spin-orbit parameters: 

£[(ss)l(33)10], £[(ss)l(33)34], £[(ss)l(33)54], £>[(ss) 1(33)56]; 

• 3 ligand-field parameters: 
£[(ss)0(33)00], £[(ss)0(33)44], £[(ss)0(33)66]. 

The generalized symmetry-adapted weak-field model with these 33 parameters is equiva- 
lent to the generalized strong- field model for f N in O [58] . 

2.6 Transition intensities 

In addition to be useful for the calculation of energy levels of a partly-filled shell ion in a 
given surrounding, the Racah's methods proved to be of considerable importance for the 
calculation of transitions between levels. We shall not develop these facets of crystal- and 
ligand-field theory here. It is enough to mention the pioneer works by Judd [85] and Ofelt 
[86] for one-photon electric dipolar transitions between split levels of the same parity (see 
also [10]). Let us also mention that the symmetry considerations developed by Bader and 
Gold [87] for two-photon electric dipolar transitions between states of opposite parities 
were reformulated in the symmetry-adapted weak- field model [59l EHl E9]. Finally, let us 
mention that irreducible tensor methods for finite groups were used for calculating the 
intensities of photoelectron spectra of partly-filled shell ion systems [90j ED [92j [93j [94]. 



3 Wigner-Racah algebra for a finite or compact group 

An important task in spectroscopy is to calculate matrix elements in order to determine 
energy spectra and transition intensities. In the case of many-fermionic systems, this can 
be done either in the Slater-Condon-Shortley approach (with determinantal states) or 
in the Dirac- Wigner-Racah approach (with states characterized by quantum numbers). 
In the Dirac- Wigner-Racah approach, one way to incorporate symmetry considerations 
connected to a chain of groups (involving symmetry groups and classification groups) 
is to use the Wigner-Racah calculus' associated with the chain under consideration. 
The Wigner-Racah calculus' or Wigner-Racah algebra' associated with a group G (or a 
chain of groups G a 3 Gr) is generally understood as the set of algebraic manipulations 
concerning the coupling and recoupling coefficients for the group G (or the head group 
G a ). This 'algebra' may be also understood as a true algebra in the mathematical sense: It 
is the (in) finite-dimensional Lie algebra spanned by the irreducible unit tensor operators 
or Wigner operators of G (or G a ) El ESI ESI E7j. We shall mainly focus here on the 
very basic aspects of the 'algebra' of the coupling and recoupling coefficients of a finite 
or compact group G. The Wigner-Racah calculus was originally developed for simply- 
reducible (i.e., ambivalent plus multiplicity-free) groups [981 EHJ Q.00J. (Let us recall that 
a group G is said to be ambivalent if each element of G and its inverse belong to a 
same conjugation class. It is said to be multiplicity- free if the Kronecker product of 
two arbitrary irreducible representations of G contains at most once each irreducible 
representation of G.) The bases of the Wigner-Racah algebra of the rotation group, a 
simply-reducible group, were introduced at the beginning of the forties by Wigner [99] 
and Racah [21 E]- In the sixties and seventies, the idea of a Wigner-Racah algebra was 
extended to an arbitrary finite or compact group [10Q 11021 1103j (see the review in [104] ) 
and started to be applied to some groups or chains of groups of interest in crystal- and 
ligand-field theory [321 H51 HO ESI ED]- Regarding molecular and solid-state physics, let 
us also mention that Koster et al. published the first complete set of tables of coupling 
coefficients for the thirty- two (single and double) crystallographic point groups |105] . Most 
of the developments concerning chains of groups were strongly influenced by a lemma due 
to Racah derived in [I] for an arbitrary chain involving finite and/or compact groups. 

We present in what follows the basic ingredients for the Wigner-Racah algebra of a 
finite or compact group in a terminology easily adaptable to nuclear, atomic, molecular, 
and condensed matter physics as well as in quantum chemistry. 

3. 1 Preliminaries 

Let us consider an arbitrary finite or compact continuous group G having the IRCs a, 
b, ■ • -. The identity IRC, often noted A or A\ or T\ in molecular physics, is denoted by 



in this section (it is noted To in Section 2). To each IRC a, we associate a unitary 
matrix representation D a . Let [a] be the dimension of D a . The a-a' matrix element of 
the representative D a (R) for the element R in G is written D a (R) aa /. (For a = 0, we have 
a = a' = 0.) The sum x a (R) — J2 a D a (R) aa stands for the character of R in D a . The 
D a (R) aa i and x a (R) satisfy orthogonality relations (e.g., the so-called great orthogonality 
theorem for D a (R) aa /) that are very familiar to the physicist and the chemist. We use 
\G\ to denote the order of G when G is a finite group or the volume J G dR of G when G 
is a compact continuous group. Furthermore, the notation j G ■ ■ ■ dR, which applies when 
G is a compact continuous group, should be understood as YIr^g • • • wnen G is a finite 
group. 

3.2 Clebsch-Gordan coefficients 

The direct product a <g> b of two IRCs a and b of G can be in general decomposed into a 
direct sum of IRCs of G. This leads to the Clebsch-Gordan series 

a®b = Q)cr(c\a®b)c (19) 

c 

where a(c\a <g> b) denotes the number of times the c IRC occurs in a <g> b. The integers 
a(c\a ® 6) may be determined through the character formula 

a( c \a ®b) = IGT 1 / xWx a (i?)x b (^ (20) 

In terms of matrix representations, (fl9~I) reads 

D a ®D b ~ 0a(c|a(g)6)L> c (21) 

c 

Therefore, there exists a unitary matrix [7 afe such that 

{U ab ) ] D a {R) ® D\R)U ab = a{c\a ® 6)D c (i2) (22) 

c 

or equivalent ly 

D a (i?) ® D 6 (i2) = a(c|a ® b)U ab D c (R){U ab ) ] (23) 

c 

for any i? in G. It is a simple exercise in linear algebra to transcribe (122"]) and (|23|) in 
matrix elements. We thus have 

E ( Ua %,, P n D a (R)aa>D b (R) m , {U ab ) aipi pldil = A(c\a ® ^S^S^R)^ (24) 

and 

DURU&Ww = J2 ( Uab ) a ^ DC ( R )^> { Ua % M ( 25 ) 

PC77' 



for any R in G. Each row index of U ab consists of two labels (a and /3) according to the 
rules of the direct product of two matrices. Similarly, two labels (c and 7) are required 
for characterizing each column index of U ab . However, when c appears several times in 
a <g> b, a third label (the multiplicity label p) is necessary besides c and 7. Hence, the 
summation over p in (125)) ranges from 1 to a(c\a <8> b). Finally in (121)1 . A(c|a <g> 6) = or 1 
according to whether as c is contained or not in a <g> 6. (Note that A(c|a <8> &) is the analog 
of A(S,L,J) used in Section 2.) 

Following the tradition in quantum mechanics, we put 

(aba(3\p C1 ) := (U ab ) a ^ (26) 

so that ( 124)) and ( T25|) can be rewritten as 

{abaplfxrfY D a (R) aa ,D b (R)^(aba' 'f3'\p V7') = A(c|a <g> b)5 pp/ 5 cc ,D c (R)^, (27) 

and 

D a (R) aa ,D b (R)pp, = ^(a6a/3|pc7) J D c ( J R) 7Y (a6a'/3'|pc7')* (28) 

PC77' 

The matrix elements (a6a/3|pc7) are termed Clebsch-Gordan coefficients (CGCs) or vector 
coupling coefficients. The present introduction clearly emphasizes that the CGCs of a 
group G are nothing but the elements of a unitary matrix which reduces the direct product 
of two irreducible matrix representations of G. As a consequence, the CGCs satisfy two 
orthonormality relations associated with the unitary property of U ab : 

^^(aba(3\pc^)*{aba(3\p'c'i) = A(c\a <g> b)5 pp >5 cc >5 y y (29) 

a/3 

and 

'y^ j (aba(3\pcy)(aba' ft' \ pay)* = Saa'bpp 1 (30) 

pc-j 

Note that ( 129)) and (130)) are conveniently recovered by specializing R to the unit element 
E of G in (127)) and (1251) . respectively. As an evident selection rule on the CCGs, it is 
clear that in order to have (aba/3\pcy) 7^ it is necessary (but not sufficient) that c be 
contained in a ® b. 

Equations (1271) and (125)) show that the CGCs are basis-dependent coefficients. In this 
regard, it is important to realize that ( 127)) and (125)) are not sufficient to define unambigu- 
ously the CGCs of the G group once its irreducible representation matrices are known. 
As a matter of fact, the relation 

(aba(3\rcy) := J ^^(abaf3\pcy)M(ab,c) pr (31) 
p 



where M(ab, c) is an arbitrary unitary matrix of dimension a(c\a®b) x a(c\a®b), defines a 
new set of CGCs since ( 12 7 jl and (1281) are satisfied by making replacements of type p — > r. 
The CGCs associated with a definite choice for the irreducible representation matrices 
of G are thus defined up to a unitary transformation, a fact that may be exploited to 
generate special symmetry properties of the CGCs. 

Various relations involving elements of irreducible representation matrices and CGCs 
can be derived from (1271) and (1281) by using the unitarity property both for the represen- 
tation matrices and the Clebsch-Gordan matrices. For instance, we obtain 

D a (R) aa/ D b (R)^(aba' 'P'\pcy') = Y^i aha P\P c l) DC { R )ii' ( 32 ) 

a'/3' 7 



J2 Da ( R )aa'(aba'^\pcj') = J2( aba ^\pn)D b (R)} f3 ,D c {R)^ Y (33) 

a' /37 

(aba'(3'\pci) = Y J ^af3\pc 1 )D\RY aal D\Ry^D c {R)^ (34) 

a/37 

for any R in G. In the situation where the elements of the irreducible representation 
matrices of G are known, Eqs. (13"2"j) . (I3"3"|) and (I3~4"|) provide us with linear equations useful 
for checking the numerical values of the CGCs of G. 

The combination of ( )28|) with the great orthogonality theorem for G yields the relation 

\G\~ l I ^(i?)^/^)^//^);^ (35) 

p 

which is useful for the calculation of the CGCs of G in terms of the elements of the 
irreducible representation matrices of G. Note that when a <g> b is multiplicity-free (i.e., 
when there is no summation on p in (135]) ). Eq. (135]) allows us to determine (aba(3\c~f) for 
all a, (3 and 7 up to arbitrary phase factors ; more precisely, we then have 

■u h s ( \c] V /2 LD a (R) aa ,D\R) m D c (R)* ,dR 

abaB 1 cV) = e lfc(ab ' c) — — — — v ' n (36) 

1 \\G\J {f G D«(R) aW D>(R)^Dc(R); y dRy/* W 



where h(ab, c) G R. 

It appears from (I32^) - (I3o]) that c does not generally play the same role as a and b 
in (aba(3\pcy). Indeed, (134"]) shows that the CGCs (afea/9|pc7) are the components of a 
third rank tensor, twice contravariant and once covariant. Therefore, (abad\pcy) does not 
generally exhibit simple symmetry properties under permutations of a, b and c. It will be 
shown in the following how the CGCs may be symmetrized thanks to a 2-aa symbol. 



3.3 The 2-aa symbol 

Let us define the 2-aa symbol through 

a b 



[a] 1/2 (6a/3a|00) (37) 

a (3 / 

The 2-aa symbol makes it possible to pass from a given irreducible matrix representation 
to its complex conjugate. This is reflected by the two relations 

(a b \ * / a 1/ \ 

a J D a (R) aa ,[ a> ^)=A(0\a®b)5 bb ,D b (R)^, (38) 



and 



a J/>(/^, Q =A(0\a®b)5 aa ,D a (R) aa , (39) 

that hold for any R in G. The proof of ( 138]) and (139]) is long ; it starts with the in- 
troduction of ( 157|) into the left-hand sides of (13"8"]) and (1391 and requires repeated use of 
relations involving the irreducible matrix representations and CGCs as well as the great 
orthogonality theorem of G. By taking R = E in ([38]) and ( 139]) . we get the useful relations 
^ / a b\* (a b'\ 

S( ,) ( a ^*«** <«» 



y — , i a b\ ( a' b \ 

J2( a)( , a) =A(0|o®6)^ aa / (41) 



and 

^ (a /?/ V«' /? 

which give back ([29]) as particular case. 

The 2-aa symbol turns out to be of relevance for handling phase problems. In this 
regard, both ([38]) and pjj lead to 



where the Frobenius-Schur coefficient 



IGr 1 / x°(i2 2 )di2 (43) 

J G 



is 1, —1, or according to as D a is orthogonal, symplectic, or complex (i.e., integer, 
half- integer or complex in Wigner's terminology). Note that 

/ b a\ (a b\ 

c« I = <U J (44) 



(3 a J \a (3 

satisfies ( 142]) . Equation ( 144]) reflects the symmetry of the matrix which enables to pass 
from the matrix D a to its complex conjugate (D a )* (cf., the Frobenius-Schur theorem). 
Thus, the 2-aa symbol plays the role of a metric tensor that transforms D a into (D a )*. It 
generalizes the Herring- Wigner metric tensor introduced for the 577(2) group (see |99j). 



3.4 The (3-aa) p symbol 

We now define the (3-aa) p symbol via 

/a b c\ . = J2[c']- 1/2 M(ba,c') plp ( C ° ' ) (baBa\p> 'H) (45) 
V« P 7/ p V7 7/ 

where M(ba, c') is an arbitrary unitary matrix. Conversely, each CGC can be developed 
in terms of (3-aa) p symbols since the inversion of ( 145]) gives 

(aba(3\pc 1 ) = [c} l / 2 J2M(ab,c); p , (\ °\ ( b °\ (46) 
p'c'y ^7 7/ VP a i / p! 

after utilization of the unitarity property of the 2-aa symbol and of the matrix M(ba, d). 

All the relations involving CGCs may be transcribed in terms of (3-aa) p symbols. For 
example, the orthonormality relations (|29|) and (|30|) are easily amenable to the form 

^-^ (a b c\ (a b c\* 

and 



.a B j J n \ a' B' 7 . 



Ef a b c\* f a b c' \ , 
„ , = A(0|a<g)&®c)<V4c'<Mc] (4? 
. V a B 7 / o V a B 7 / . 



af3 

Along the same line, the introduction of ( 1475]) into (128]) yields 



, ^-^ f a b c\ fa b c \ * 

D«(R) aa ,D b (R)w=J2lc}( a ) D%R)^,[ , (49) 

^ V" P 7/p Va /3 7 /„ 

which in turn leads to 

^-^ f a b c\* , f a b c' \ 

= A(0|o (8) 6 (8) ^^^[c]- 1 ^^^)^, (50) 

owing to the orthogonality relation f )48|) . Equations ( I49p and f )50|) hold for any element i? 
in G. As a check, note that for R = E, they can be specialized to ( |47|) and (]48p . 

Relation (149]) and its dual relation ( 150]) show that -D a , -D 6 and D c present the same 
variance. This may be precised by 

(a b c \ ^-^ / a b c\ , 

U ? y), = £U , J^W-^W («) 

which shows that the behavior of the (3-aa) p symbol under permutations of a, 6 and c 
should be easier to describe than the one of the CGC (abaB\pcy). This is reflected by the 
following relation (to be compared to (135]) ) 



1 f $. ( a b c\ ( a b c \ * 

G - 1 / D a (R) aa ,D b (R)^D c (R) 7Y dR = J2( a )(/*/) ( 52 ) 
/g „ \a /3 jj \a' B 7/ 



which may be proved directly by combining f H9|) with the great orthogonality theorem 
for the G group. When the triple direct product a Cg> 5 £g> c contains the identity IRC of G 
only once (i.e., when there is no label p and no summation in ( 152]) ). Eq. ( 1521) shows that 
the square modulus of the 3-aa symbol is invariant under permutation of its columns. 
In this case, we may take advantage of the arbitrariness of the matrix M in (l3Tj) or 
( 145]) to produce convenient symmetry properties of the 3-aa symbol under permutations 
of its columns. By way of illustration, let us mention the following result [99]: For G 
simply reducible, it is possible to arrange that the numerical value of the 3-aa symbol be 
multiplied by the phase factor (—l) a+b+c , with (— l) 2x = c x , under an odd permutation of 
its columns ; consequently, the numerical value of the 3-aa symbol remains unchanged 
under an even permutation of its columns (since c a c b c c = 1). 

To close this subsection, we note that the (3-aa) p symbol constitutes a generalization 
to the case of an arbitrary finite or compact group of the 3-jm symbol introduced by 
Wigner for simply reducible groups (in particular for the rotation group) [99] and of the 
V symbol introduced by Fano and Racah for the SU(2) group [7] (the V symbol is a 
symmetrized version of the V symbol defined by Racah [2]). 

3.5 Recoupling coefficients 

We now define two new coefficients: 

(a(bc)p bc c bc p'd'5'\(ab)p ab c ab cpd5) := ^ (aba(3\p ab c ab j ab )(c ab c-f ab -f\pd5) 

x {bc(3^\p bc c bc ^ bc )*(ac bc a^ bc \p'd'5')* (53) 

and 

{{ac)p ac c ac (bd)p bd c bd p'e'e'\(ab)p ab c ab (cd)p cd c cd pee) 
■ = ^ ^2 ( aba (3\PabC a blab)(cd75\p cd c cd j cd )(c ab c cd j ab 'y cd \pee) (54) 

-Jablcd laclbd 

x (aca^\p ac c ac 'y ac y(bdf3d\p bd c bd -f bd y(cacC bd 'yaclbd\p f e'e'y 

The introduction in these definitions of and the use of the great orthogonality theorem 
for G leads to the properties 

(a(bc)p bc c bc p'd'5' \ (ab)p ab c ab cpdS) 
= 8 dd/ 5 S 6> [d]' 1 y^ j (a(bc)p bc c bc p'd5\(ab)p ab c ab cpd5) (55) 

and 

{{ac)p ac c ac (bd)p bd c bd pe'e'\(ab)p ab c ab (cd)p cd c cd pee) 

= bee'bee'W 1 ^((ac) p ac c ac {bd) p bd c bd p' ' ee\{ab) p ab c ab {cd) p cd c cd pee) (56) 



so that the recoupling coefficients defined by (1531 and fl54l) are basis-independent (i.e., 
they do not depend on the labels of type a) in contrast with the coupling coefficients 
(abaf3\pc r y). 

By using the orthonormality of the CGCs, it can be shown that the CCGs occurring 
in Eqs. (|53|) and (1541) can be moved from the right hand side to the left hand side in such 
a way to produce new relations for which the total number of CGCs remains equal to 4 
and 6, respectively. Repeated actions of this type lead to orthonormality relations for the 
recoupling coefficients (l5"3"j) and (15"4"|) . 

In a way paralleling the passage from the coupling coefficients to the (3-aa) p symbol, 
one can define (6-a)4 P and (9-a)s p symbols from the recoupling coefficients defined by 
fl5I])-(ES]). The defining expressions (6-a) 4p and (9-a) 6p symbols are very complicated and 
not especially instructive in the case of an arbitrary compact group G. Hence, they shall 
be omitted as well as the defining expressions for higher (3N-a)2N P symbols corresponding 
to the recoupling of N > 4 IRCs. Finally, note that the recoupling coefficients and their 
associated (3N-a) 2 N P symbols, N > 1, for a G group can be connected to other basis- 
independent quantities, viz., the characters of G |1U1[ 1106] . 

3.6 Irreducible tensorial sets 

Let {|Taa) : a = 1, 2, . . . , [a]} be a basis for the irreducible matrix representation D a of 
G. The vectors \raa) are defined on a unitary or pre-Hilbert space £ (indeed, a Hilbert 
space in the quantum-mechanical applications) and there exists an application R \-¥ Pr 
such that 



for any R in G. Following the work by Fano and Racah [7J on the SU(2) group, we refer the 
set {|raa) : a = 1, 2, . . . , [a]} to as an irreducible tensorial set (ITS) of vectors associated 
with D a . The label r may serve to distinguish different ITSs of vectors associated with the 
same irreducible matrix representation D a . (In practical applications, this label consists 
of various quantum numbers arising from nuclear, or atomic or molecular configurations.) 
In this connection, note the following standardization: It is always possible to arrange 
that { | Tact) : a = 1, 2, . . . , [a]} and {|r'aa) : a = 1, 2, . . . , [a]} span the same matrix 
representation D a rather than two equivalent representations. We shall assume that such 
a standardization is always satisfied. 

From two ITSs {\r a aa) : a = 1, 2, . . . , [a]} and {\ubf3) : (3 = 1, 2, ... , [6]}, we can 
construct another ITS of vectors. Let us define 



[a] 




(57) 



Q' 



r'=l 




(58) 



Then, as a simple corollary of ([28]) . the set {\T a Tt,abpcy) : 7 = 1, 2, . . . , [c]} can be shown 
to be an ITS associated with D c . 

In a similar way, let us consider a set {T£ : a = 1,2, . . . , [a]} of (linear) operators 
defined on £ and such that 

[a] 

PrT^p- 1 = J2 T Z' Da ( R )*'« ( 59 ) 

a'=l 

for any R in G. This set is called an ITS of operators associated with D a . We also say that 
this set defines an irreducible tensor operator T a associated with D a . Note the implicit 
standardization: The sets {T^ : a = 1, 2, . . . , [a]} and {U^ : a = 1, 2, . . . , [a]} span the 
same matrix representation D a rather than two equivalent representations. 
In full analogy with flo^]) . we define 

{X a 3 v *yc ._ £ T Q "^( a6a/3 | pc7 ) ( 6 o) 

a/3 

from the two ITSs {T^ : a = 1, 2, . . . , [a]} and {U b : j3 — 1, 2, . . . , [&]}. As a result, the set 
{{T a £g> U 6 }^ c : 7 = 1, 2, . . . , [c]} is an ITS of operators associated with .D c . We say that 
{T a ® U 6 } is the direct product of the irreducible tensor operators T a and XJ b . Observe 
that this direct product defines a tensor operator which is reducible in general. Equation 
(I6"0"|) gives the various irreducible components of {T a ® XJ b }. 

3. 7 The Wigner-E chart theorem 

The connection between most of the quantities introduced up to now appears in the 
calculation of the matrix element (T'a'a'\T^\Taa), the scalar product on S of the T^\raa) 
vector by the {r'a'a') vector. By developing the identity 

(r'a'a'\T b \Taa) = (r'a'a'\P R P R T b P R l P R \raa) (61) 

we get, after some manipulations, the following basic theorem. 

Theorem 1 (Wigner-Eckart's theorem). The scalar product (t' a' a'\Tp\Taa) can be 
decomposed as 

(r'a'a'\T^\Taa) = yj(r'a'| \T b \ \ra) p 

p a" a" 

Alternatively, (162]) can be cast into the form 

(T'a'a'\T b \Taa) = [a'}'^ ^(r'a'| \T b \ \ra) p (aba(3\pa'a'Y (63) 

p 

with 

{T'a'\\T b \\Ta} p := ^ M(ab, a') Pp ,(T'a'\ \T b \ |ra)^ (64) 
p' 




where M(ab,a') is an arbitrary unitary matrix (cf., (|45j) and (1461) ). 

In the summation-factorization afforded by f )62|) or (163]) . there are two types of terms, 
namely, the (3-aa) p symbols or the CGCs (aba (3\ pa' a') that depend on the G group only 
and the so-called reduced matrix elements {r'a'\ \T b \ \ra) p or (r'a'\ \T b \ \ra) p that depend 
both on G and on the physics of the problem under consideration. The reduced matrix 
elements do not depend on the 'magnetic quantum numbers' (a', (3 and a) and therefore, 
like the recoupling coefficients, are basis- independent. We then understand the interest 
of the recoupling coefficients in applications: The reduced matrix elements for a com- 
posed system may be developed as functions of reduced matrix elements for elementary 
systems and recoupling coefficients. In this direction, it can be verified that the matrix 
element (r^a'b 1 p'c ,r y'\{T d ®\] e } a ^\r a T b abpc~i) can be expressed in terms of the recoupling 
coefficients defined by (154"]) and (156"]) . 

Equations (162]) and (163]) generalize the Wigner-Eckart theorem originally derived by 
Eckart for vector operators of the rotation group |107j . by Wigner for tensor operators 
of the rotation group |1U8] and of simply reducible groups [22], an d by Racah for tensor 
operators of the rotation group [2]. 

A useful selection rule on the matrix element (r'a'a'\Tl\Taa) immediately follows from 
the CGCs in fl63|) . The latter matrix element vanishes if the direct product a® b does not 
contains a'. Consequently, in order to have (r'a'a'\T^\raa) ^ 0, it is necessary (but not 
sufficient in general) that the IRC a' be contained in a <g) b. 

As an interesting particular case, let us consider the situation where b is the identity 
IRC of G. This means that the operator H = Tq is invariant under G (see (159|) ). Equation 
( 16"3"]) can be particularized to 

(r' 'a' 'a'\H\raa) = 5 aa '5 aa '(r'a\\T \\Ta) (65) 

where the index p is not necessary since a <E> = a. The Kronecker deltas in ( 165]) show 
that there are no a' -a and/or a! -a. mixing. We say that a and a are 'good quantum 
numbers' for H . The initial and final states have the same quantum numbers as far as 
these numbers are associated with the invariance group G. The invariant H does not mix 
state vectors belonging to different irreducible representations a and a'. Furthermore, it 
does not mix state vectors belonging to the same irreducible representation a but having 
different labels a and a'. 

It is very important to realize that phase factors of type (— l) a , (— l) a ~ Q and (— l) a+6+c 
do not appear in f )62|) and fl63|) . Indeed, the present exposure is entirely free of such phase 
factors, in contrast with other presentations. As a matter of fact, in many works the 
passage from the Clebsh-Gordan or unsymmetrical form to the (3-aa) p or symmetrical 
form of the coupling coefficients involves unpleasant questions of phase. This is not the 
case in (T4"5]) and (146]) . Such a fact does not mean that (14"5]) and (146]) as well as other 



general relations are free of arbitrary phase factors. In fact, all the phase factors are 
implicitly contained in the matrices M, the 2-aa symbols and the (basis- independent) 
Frobenius-Schur coefficient. 

3.8 The Racah lemma 

We have already emphasized the interest of considering chains of groups rather than 
isolated groups. Let us now denote G as G a and let Gr be a subgroup of G a . In this case, 
the labels of type a, that occur in what precedes, may be replaced by triplets of type aY^. 
The label of type Y stands for an IRC of the group Gr, the label of type 7 is absolutely 
necessary when [r] > 1 and the new label of type a is a branching multiplicity label to 
be used when the V IRC of Gr is contained several times in the a IRC of the G a head 
group. (The 7 label is an internal multiplicity label for Gr and the a label is an external 
multiplicity label inherent to the restriction G a — > Gr-) Then, the (ciia 2 aia 2 |paa) CGC 
for the G a group is replaced by the (aia 2 aiFi7ia 2 r 2 7 2 |paar7) CGC for the G a group in 
a G a D Gr basis. We can prove the following theorem. 

Theorem 2 (Racah's lemma). The CGCs of the G a group in a G a D Gr basis can be 
developed according to 

(aia 2 «iri7 1 a 2 r 2 72|paar7) = ^^^27172 ^7) (ai 0:^1 + a 2 a 2 T 2 \paaT)p (66) 

P 

where the (Fir 2 7i7 2 |/3r7) coefficients are CGCs for the Gr group considered as an isolated 
group and the (aiaiTi + a 2 a 2 r 2 |paar) ( g coefficients do not depend on ji, 7 2 and 7. 

The proof of Racah's lemma was originally obtained from Schur's lemma [I]. However, 
the analogy between f )62|) . fl63|) and f )66|) should be noted. Hence, the Racah lemma for 
a G a D Gr chain may be derived from the Wigner-Eckart theorem, for the G a group 
in a G a D Gr basis, applied to the Wigner operator, i.e., the operator whose matrix 
elements are the CGCs. The (aiO^Tx + a 2 a 2 r 2 |paar) (9 in the development given by 
( 166]) are sometimes named isoscalar factors, a terminology that comes from the SU (3) D 
U(l) <g> SU{2) chain used in the eightfold way model of subatomic physics. 

From a purely group-theoretical point of view, it is worth to note that Racah's lemma 
enables us to calculate the CGCs of the Gr subgroup of G a when those of the G a group 
are known (see for example |109] and references therein). In particular, for those triplets 
(rir 2 r) for which Ti £g> T 2 contains V only once, the CGCs (FiF 2 7i7 2 |r7) are given by a 
simple formula in terms of the CGCs of G a - 

The summation-factorization in f[6"6"j) can be applied to each CGC entering the defini- 
tion of any recoupling coefficient for the G a group. Therefore, the recoupling coefficients 
for G a can be developed in terms of the recoupling coefficients for its subgroup Gr [2U ED] . 



3.9 Illustrative examples 

3.9.1 The 5*7(2) group in a SU(2) D U(l) basis 

As a first example, we take G a = ££7(2) and Gt = U(l) where SU(2 and £7(1) are 
the universal covering groups or, in the terminology of molecular physics, the 'double' 
groups of the proper rotation groups R(3) ~ SO (3) and 72(2) ~ 5*0(2), respectively. 
In this case, a = (?) where j is either an integer (for vector representations) or a half- 
of-an-odd integer (for spinor representations), aT^j = m ranges from — j to j by unit 
step, and D a (R) aa i can be identified to the element D^\R) mm i of the well-known Wigner 
rotation matrix of dimension [j] = 27 + 1. The matrix representation D^' corresponds 
to the standard basis {\j, m) : m — j, j — 1, . . . , — j} where \j, m) denotes an eigenvector 
of the (generalized) angular momentum operators 7 2 and J z . (For j integer, the label 
I often replaces j.) The labels of type m clearly refer to IRCs of the rotation group 
Coo ~ -R(2). Therefore, the basis {\j, m) : m = j, j — 1, . . . , — j} is called a R(3) D R(2) or 
SU{2) D Z7(l) basis. Furthermore, the multiplicity label p is not necessary since SU{2) 
is multiplicity- free. Consequently, the (real) CGCs of 577(2) in a SU{2) D U(l) basis are 
written (ji^^i^ljTn)- They are also called Wigner coefficients. 

In view of the ambivalent nature of SU(2), the 2-aa symbol reduces here to 

" ' ' <\,r ( 3 ,) (67) 

\ 777. 777, J \m 777, 

We can take 



777, 777,' 



(-l) J+m 5(m', -m) (68) 



where (— l) jf+m 5(m', — m) is a component of the 1—jm Herring- Wigner metric tensor (in 
the Edmonds normalization |110j ). Then, the introduction of (l67|l and (|68|) into ( 145|) 
for the SU(2) D t7(l) chain shows that the 3-aa symbol identifies to the 3-jm Wigner 
symbol 

( 31 32 33 ) ■= (2j 3 + l)-H-iy 3 - m3 - 232 (j2jim 2mi \j 3 , -m 3 ) (69) 
\mi ni2 rns J 

provided we chose M(j 2 ji, js) = (— l) 2jfl . Such a choice ensures that the 3-jm symbol is 
highly symmetrical under permutation of its columns. 

In the £77(2) case, the (6-a) 4p and (9-a)Q P symbols may be chosen to coincide with the 
6-7 Wigner (or W Fano-Racah) symbol and the 9— j Wigner (or X Fano-Racah) symbol, 
respectively. More precisely, we have 

/; 323 ' ) ■= (-l) il+i2+i3+ i(2ii 2 + l)(2i 23 + l)]-^ 

J3 312 32 J 

x {h{hh)hz3^\{hh)3i2h3m) (70) 




and 



:= {(2 Jl2 + l)(2j 3 4 + l)(2j'i3 + l)(2j 24 + I)]"' 

x {{jih)3iz{323i)32ijm\ (ji^j^OW^W'rn) (71) 

in terms of recoupling coefficients (cf., (155]) and (!56|) ). 

Finally, for a = (k), the T a ITS coincides with the T^ k ' irreducible tensor operator 
of rank k (and having 2k + 1 components) introduced by Racah. We denote by Tg k ^ the 
components of T( fc ) in a SU{2) D ?7(1) basis. 

All the relations of subsections 3.1-3.7 may be rewritten as familiar relations of angular 
momentum theory owing to the just described correspondence rules. For example, (1381) 
or f[3"9"j) and f[6"2"l) can be specialized to 

D U) (Ry mm , = (-l) m - m 'D«\R)_ m> _ m , (72) 

and 

(rYm'lTf \rjm) = (-lY'- m ' ( J \ * J Vx'/IIT^I^) (73) 

\ — m q m J 

respectively. For more details, the reader should consult the textbooks in Refs. [7J II lUj 
(see also |g[iniin]). 



3.9.2 The SU{2) group in a SU{2) d G* basis 

We now consider the case G a = SU(2) and Gr = G*, where G* is isomorphic to the 
double group of a point (proper) rotation group G. Then, we have a = (J) and we take 
aT r ) = aT r y for the labels a and aY^ of Section 3.8. This will be clarified below. 
1 - The restriction of SU{2) to G* 

Each IRC (j) of SU{2) can be decomposed into a direct sum of IRC's of G*\ 

(i) = I>(r|j)r (74) 

r 

where 

a(T\j) = IGT 1 [ dRx T {RTx {3 \R) (75) 
Jg* 

stands for the multiplicity of the V IRC of G* in (j). In terms of unitary matrix repre- 
sentations, this means that 



DO ~0a(r|j)D r 
r 



(76) 



In other words, there exists a unitary matrix [7 J such that 

(UyD®(R)U* = Q)a(T\j)D r (R) (77) 
r 

holds for any R in G*. This leads to 

Y,{Mj^irD^\R) mml {jm'\ja'T'i) = S^.D 7 {R)^, (78) 

mm' 

or 

D^\R) mm , = (jmljaT^D^R^ijm'ljaTj')* (79) 

aT'fy' 

for any R in G*. In (ITS]) and (IT9]) . (jm\jaT , y) denotes an element of the matrix W: 

(jm\jaTj) := (W) m<ari (80) 

The label a (cf., the column index aT^y of £P) is a branching multiplicity label indispens- 
able when T appears more than once in (j). Note that the unitary property of the matrix 
W corresponds to R = E, the unit element of G*, in (ITS]) and (]T9l : 



^2(jm\jaT>yy(jm\jaT'>y') = cW<5rT"V/ (81) 

m 

or inversely 

mm' V° I 

Observe that (ITS]) and (IT9]) are note sufficient for determining the reduction coefficients 
{jm\jaT r )) once the irreducible representation matrices of G* and SU(2) are known since 
the coefficients 

(jm\jbTj) := J2(Mj^l)M ab (83) 

a 

where M is an arbitrary unitary matrix satisfy f ]TS]) and f ]T9"]) with the replacement a — > b. 
Nevertheless, ( ITS]) and (IT9]) lead to systems that may be useful for the calculation of the 
{jm\jaT r )) coefficients. 

2 - Irreducible tensorial sets 

From the ITS of vectors {\Tjm) : m = j,j — 1, . . . , — j} associated with we define 

\rjaTj) := ^ \ T j m )U m \j aT l) ( 84 ) 
in 

Equation (T7U|) allows us to show 

P R \rjaYi) = \rjaTi)D r (R)^ (85) 



for any R in G*. Similarly, from the ITS of operators {Tg k ^ : q = k,k — 1, . . . , — k} 
associated with D(% we define 

Ti kaT) = : = 53T«(*g|*or 7 ) (86) 

so that 

P B T^ T) P R l = J2 T^ aV) D T {R)^ (87) 

7' 

holds for any R in G* . 

At this point, it is important to remark that f )84p and (1861) provide us with ITSs both 
for SU(2) and G*. Indeed {\rjaT^) : 7 ranging } is an ITS of vectors spanning the matrix 
representation D r of G* while {\TjaT / -f) : aTj ranging } is an ITS of vectors spanning the 
matrix representation T>^> of 577(2) defined by 

T> ij \R)ar 7 ,a'r>y := Yj(Mja£lY^\R)rn m <{jrn'\ja'T'i) (88) 

mm' 

for any R in SU(2). A similar remark applies to the sets {T^ kar ^ : 7 ranging } and 
{ T ir 7 : aF 7 ranging}. 

3 - Wigner-Eckart theorems 

As an important consequence of the latter two remarks, we may apply the Wigner- 
Eckart theorem either for the group .577(2) in a SU (2) D G* basis or for the group G* in 
a G* C S77(2) basis. For G* in a G* C SU(2) basis, (J62]) gives 

(rijia 1 r 1 7 1 |Tf ar )|r 2j2 a 2 r 2 7 2 ) = ^(r 1 j 1 a 1 r 1 ||T( fear )||r 2 j 2 a 2 r 2 ) p 



x 

1. v.; 



E( r; Tl )( T T > r; ) («») 

v 7i 7i / V 7 72 7i 



For SC/(2) in a SC/(2) D G* basis, we can combine fH, ® and QBE) to obtain the 
compact formula 

(r 1 J 1 a 1 r 1 7i|T a { ^|r 2 j 2 a 2 r 2 7 2 ) = (r 1 j 1 ||T( fc )||r 2 j 2 )/f * ^ M (90) 

\airi7i a 2 r 2 7 2 6^7/ 

where the / symbol is defined by 

/ f * t 2 t) ■= (-l) 2 "(2ji + l)- 1/2 (j2fca 2 r 2 7 2 ar 7 |j 1 a 1 r 1 7 1 )* (91) 

\ ail 171 a 2 l 2 7 2 al7/ 

in function of the CGC 

(jij2aiFi7ia 2 r 2 7 2 |jar7) := ^ (ji m i|jiai r i7i)* 

x (j2m2|j2a 2 r 2 7 2 )*(jij 2 mim 2 |jm)(jm|jar7) (92) 




of SU(2) in a SU(2) D G* basis [26]. 



There are many advantages to use (190]) rather than fl89|) . In fl89|) . both the reduced 
matrix elements and the coupling coefficients (cf., the 2-aa and (3-aa) p symbols) depend 
of the symmetry group G* . Furthermore, the factorization offered by (1891) is not complete 
in view of the summation over the multiplicity label p. On the other side, the matrix 
element given by (1901) factorizes in two parts: a coupling coefficient (cf., the / symbol) for 
the SU(2) D G* chain and a reduced matrix element which does not depend of the group 
G*. This maximal factorization takes place even in the case where G* is not multiplicity- 
free. The reduced matrix elements in (l90j) applied to complex systems either are obtainable 
from tables or can be calculated from Racah's method in terms of recoupling coefficients of 
SU(2), coefficients of fractional parentage, and elementary reduced matrix elements. The 
main calculation to do when dealing with (|9T)|) most of the time concerns the / geometrical 
coefficient, a quantity which is independent of the additional quantum numbers t\ and 
T2 and which remains invariant when the tensor operator T k is replaced by any tensor 
operator U fc of the same rank. 

The calculation of the / coefficients defined by (|9~Tj) and (l9"2l touches a simple problem 
of symmetry adaptation. In fact, the determination of the symmetry-adapted CGCs f )92|) 
require the knowledge of the reduction coefficients flHOj) . These reduction coefficients 
are the expansion coefficients of symmetry adapted functions (cf. (184"j) ) or symmetry- 
adapted operators (cf. (I56"j) so that their calculation may be achieved by numerous means 
(resolution of linear systems like (175]) or ( I79|) . projection operator techniques, . . . ). 
4 - The / symbol 

Equation ( 191 j) shows that the behavior of the / symbol under the interchange of its 
first and second columns is not easy to describe. The / symbol may be symmetrized 
owing to the introduction of the 1-jaTj symbol 



where the l—jm symbol is defined by (I6T|) and (155|) . The / or 3-jar7 symbol defined 
through 



VaJY}^ a 2 r 2 72 ailYyi/ 
then exhibits a high (permutational) symmetry since a simple developement of (|94|) leads 
to 




(93) 




(94) 




(95) 



(see [26] for the properties - selection rules, orthogonality relations, etc. - of the / and / 
symbols) . 

For G* = U(l) the / symbol and the l—jaT'y symbol reduce to the 3-jm Wigner 
symbol and to the 1—jm Herring- Wigner symbol, respectively. The l—jaT^ and / symbols 
are thus 2-aa and 3-aa symbols as defined in Sec. II (with a — > j and a — > GT7), 
respectively, for the group SU{2) in a SU{2) D G* basis. The properties (existence 
conditions, selection rules, symmetry properties, orthogonality properties, . . . ) of the / 
(and /) symbols can be deduced from the ones of the 3-jm symbols and the JJi matrices 
and have been discussed at length elsewhere [26] [271 128] . Let us simply mention that, 
by applying Racah's lemma, the / symbol can be developed as a linear combination of 
(3-T7)p according to 

/( ■" n h ) >-./(( h h h )) ( Ti T ° T:i ) (96) 

\cnri7i a 2 r 2 7 2 a 3 r37 3 / Wairi a 2 r 2 a3r 3 // p \7i 72 73 J p 

where the /((...)) reduced coefficient is independent of 71, 7 2 et 73. 

We are now in a position to enunciate correspondence rules for passing from the 
Wigner-Racah algebra of SU(2) in a SU(2) D U(l) basis (i.e., in the {jm} scheme) to 
the Wigner-Racah algebra of SU{2) in a SU{2) D G* basis (i.e., in the {jaT^} scheme): 
All the m- or g-dependent quantities are replaced by the corresponding al^-dependent 
quantities while the basis-independent quantitities (like 6-j and 9-j symbols) are un- 
changed. More precisely, we have 



DV\R) mm . 


— >■ 




Uij2m 1 m 2 \jm) 


— >• 


(jij2Ctiri7ia 2 r 2 7 2 |jar7) 


1—jm symbol 


— >• 


1— aT'-f symbol 


3— jm symbol 


-> 


/ symbol 


3(n — 1)— j symbol 


->• 


3(n — 1)— j symbol 


rjm) 


— >■ 


\rjaV-f) 


T (k) 
1 


— >■ 


T (k) 



(97) 



(see [2Z] for more details). 

3.9.3 The G* group in a G* c SU{2) basis 

1 - The general case 

Equations fl32|) - fl36|) were used in numerous works for calculating coupling coefficients 
and V symbols of subgroups of SU{2). (Following Griffith [15], the (3-aa) p symbols 
of a group of molecular interest are referred to as V symbols in what follows.) We now 
describe an alternative method for calculating the V coefficients of a subgroup G* of SU (2) 



as renormalized / coefficients of the £77(2) D G* chain. This method combines three basic 
ingredients scattered in various (implicit or explicit) approaches starting with the pioneer 
works by Tanabe, Sugano and Kamimura: the concept of quasi angular momentum, the 
definition of the / symbol and renormalization techniques. For the purpose of simplicity, 
we shall limit ourselves to a multiplicity-free group G* but it should be noted that the 
method may be extended to an arbitrary subgroup of SU{2). 

Given the V IRC of G* , let (j(T)) or simply (j) the IRC of SU{2) that contains V once 
and only once. Thus, j is the smallest value of j for which er(r|?) = 1. The value j refers 
to a quasi angular momentum (see |109] too). In the multiplicity-free case where the 
identity IRC of G* appears once and only once in the triple direct product ri <S> T 2 <8> T 3 , 
there is no need for the internal multiplicity label p in the 3-r7 or V symbol. Therefore, 
let us put 



V 



7i 



r 2 

72 



r 3 

73 



x(rir 2 r 3 )/ 



E 



.717273 



31 32 33 

Ti7i r 2 7 2 r 3 7 3 J 

3i 32 33 

Tiii r 2 7 2 r 3 7 3 



-1/2 



(9? 



where x(rir 2 r 3 ) is an a arbitrary phase factor that depends on Ti, T 2 and T 3 only. It 
can be verified by repeated application of (19T)]) that the V symbol defined by (I98p satisfies 
( I49p and f loU]) for G*. Consequently, the V symbol is nothing but a 3-r7 symbol for the 
G* group compatible with the choice implicitly assumed through ( 195]) with j = j(T) for 
the representation matrices D r . 

For the sake of simplifying calculations with (198 p . it should be noted that 



E 



717273 



3i 32 33 

Fiji r 2 7 2 r 3 7 3 



w E 

all ~/ k except % 



f 



3i 32 33 

Txji r 2 7 2 r 3 7 3 



(99) 



for % = 1, 2 or 3. In addition, if two of the three T's are equivalent to two of the 
corresponding three (j)'s, the right-hand side of (199]) can be simplified and (198]) takes a 
simple form. For instance, in the case (ji) = Ti and (j 2 ) = T 2 , (198]) becomes 

ri r 2 r 3 ' 



v 



x(r 1 r 2 r 3 )[r 3 ]- 1 / 2 (2j 3 + i) 1 /V 



31 32 33 
Tiii T212 r 3 7 3 



(100) 



. 7i 72 73 . 
which is very simple to handle. 

The main advantages of the method based on ( 19"%]) - 01 00]) for calculating the V coeffi- 
cients of G* may be seen to be the following. First, the calculation is easy in the sense 
that the V coefficients are deduced from a minimal set of / coefficients which are readily 
calculated (by hand or with the help of a computer) from (195]) . The thus obtained V 
coefficients of the G* group are simple linear combinations of 3-jm coefficients for the 
SU{2) D U(l) chain. Second, such a method allows us to work with bases of interest for 



molecular physics and quantum chemistry. In this respect, we may use in (1951) reduction 
coefficients (jm\jaT'j) corresponding to Cartesian p, d and / spin-orbitals or correspond- 
ing to a chain of groups (for instance, the SU(2) D O* D D\ D D\ tetragonal chain or 
the SU(2) D O* D D% D C| trigonal chain). Third, it is possible to transfer some of the 
features (formulas, symmetry properties, . . . ) of the 3-jm symbol from the SU(2) D U(l) 
standard chain to the V symbol of G*. For example, the permutation symmetry properties 
of the V symbol can be chosen to be essentially the ones of the S—jm symbol. In fact, by 
choosing a^ril^Ts) invariant under the 3! permutations of its arguments, the V symbol 
given by fl98l) - (U00p is multiplied by (— l)-? 1 ( ri )+-? 2 ( r2 )+i 3 ( r3 ) under an odd permutation of 
its columns so that is is invariant under an even permutation. 
2 - Application to the octahedral group 

As an illustration, we consider the case where G* is the O* double octahedral group 
and limit ourselves to the determination of the V coefficients of the O octahedral group. 
Therefore, we can replace 577(2) D O* by 5*0(3) D O. The restriction of 5*0(3) to O 
yields 



jXA 1 ) = 0, j(A 2 ) 



(101) 



where Ai, A 2 , E, T\ and T 2 denote the various IRCs of O. In view of the permutation 
symmetry properties of the V symbol, there are a priori 39 independent V coefficients 
to be calculated for the O group. The \]T r y) vectors (the label a is not necessary here) 
required for calculating these coefficients are given by 



OAlOi) 
3A 2 a 2 ) 

2E6) = 

lTtx) = 



[0,0) 



3,2)-|3,-2)] 



lT x z) 
2T 2 x) 

my) 

2T 2 z) 



1 



|2,0), |2Ee) = — [|2,2) + |2,-2)] 



~V2 



ii, r 



my) = —[[i,i) + (i 



-i)] 



(102) 



z|l,0) 



_[|2,1) + |2,-1)] 
^[|2,1) -12,-1)] 
-^[|2,2)-|2,-2)] 



in terms of spherical basis vectors \j,m) (the generic symbol 7 is a\ for A\\ a 2 for A 2 ; 
6 and e for E\ x, y and z for T\\ and x, y and z for T 2 ). The 39 independent V 
coefficients are then easily calculated from ( 195]) and (]98]) - (11021) . They are of course all 



real if we replace the pure imaginary number i by 1 in f |102j) . In the case i = a/— 1, it is 
possible to decrease the number of independent V coefficients by conveniently choosing 
the x(r!r 2 r3) phase factors. Along this line, by taking i = y—T and x(r 1 r 2 r 3 ) = 1 
except x(ET 2 T 2 ) = x(TiT]Ti) = x(TiTiT 2 ) = x{T 2 T 2 T 2 ) = —1, the reader will verify that 
Eqs. ( 195]) and ( l98]) - (ll02p lead to the real numerical values obtained by Griffith [45] for the 
V coefficients of O in his real tetragonal component system. 

It should be noted that each V coefficient calculated from (1951) and (l9~5]) - (|102p can be 
reduced (up to a multiplicative factor) to a single 3-jm coefficient for the 5*0(3) D SO (2) 
chain. We thus foresee that some properties of certain 3-jm symbols for the SU(2) D U(l) 
chain may be derived by looking at some properties induced by a subgroup of SU(2). As 
an example, we have 

(A 2 A 2 E\ _/ 3 3 2 \ / 3 3 2\ 
V \a 2 a 2 e)~ f \A 2 a 2 A 2 a 2 E6 ) ~ ~ \-2 2 Oj (iUo) 
It is clear that the value of the V coefficient in fl 1 3 [) is zero since the A 2 <g> A 2 <g> E triple 
Kronecker product does not contain the A\ IRC of the O group. As a consequence, the 

3- jm symbol in (I103p corresponding to the SU(2) D U(l) chain vanishes (owing to a 
selection rule for O) in spite of the fact that the (trivial and Regge) symmetry properties 
for SU{2) D U(l) do not impose such a result. 

To close Section 3, it is to be mentioned that diagrammatic methods initially developed 
for simplifying calculations within the Wigner-Racah algebra of the rotation group |112j 
where extended to the case of a finite or compact group [1131 IH4[ 1115] . Note also that 
considerable attention was paid in the nineties to the Wigner-Racah calculus for a q- 
deformed finite or compact group (see |116j for some general considerations on this subject 
and [1171 1118], 1119] for some developments on U q (su(2)) and U q (su(3))). 

4 Contact with quantum information 

4- 1 Computational basis and standard SU(2) basis 

In quantum information, we use qubits which are nothing but state vectors in the Hilbert 
space C 2 . The more general qubit 

:=co|0) + ci|l), c eC, dGC, |c | 2 + | Cl | 2 = 1 (104) 

is a linear combination of the vectors |0) and |1) which constitute an orthonormal basis 

5 2 :={|0),|1)} (105) 

of C 2 . The two vectors |0) and |1) can be considered as the basis vectors for the funda- 
mental IRC of SU(2), in the SU{2) D U(l) scheme, corresponding to j = 1/2 with 

|0) = |l/2,l/2), |1> = |l/2,-l/2) (106) 



More generally, in dimension d we use qudits of the form 

d-l d-1 



\ipd) •= y^c n |n), c n e C, n = 0, 1, . . . ,d - 1, ^|c n | 2 =l (107) 



n=0 n=0 

in terms of the orthonormal basis 



B d :={\n) 171 = 0,1,..., d-l} (108) 



of C d . By introducing 



m:=n _i(rf_i) ; |j, m ):=|d-l-n) (109) 

the vectors \n) can be viewed as the basis vectors for the (j) IRC of SU (2) in the SU (2) D 
U(l) scheme. In this scheme, the \j,m) vector is a common eigenvector of the Casimir 
operator J 2 (the square of an angular momentum) and of a Cartan generator J z (the z 
component of the angular momentum) of the su(2) Lie algebra. More precisely, we have 
the relations 

J 2 \j,m) = j(j + l)\j,m), J z \j,m) = m\j,m) (110) 

which are familiar in angular momentum theory. In other words, the basis Bj, known in 
quantum information as the computational basis, can be visualized as the SU(2) D U(l) 
standard basis or angular momentum basis 

B 2 j+i ■= {\j,m) :m = j,j-l,..., -j} (111) 

with the correspondence 

\0) = \j,j), |1> = U,i - 1>, \d-l) = \j,-j) (112) 

between qudits and angular momentum states. 

4-2 Nonstandard SU(2) basis 

We are now in a position to introduce nonstandard SU(2) bases which shall be connected 
in the next subsection to the so-called mutually unbiased bases (MUBs) of quantum 
information. As far as the representation theory of SU(2) is concerned, we can replace 
the set {J 2 , J z } by another complete set of two commuting operators. Following [120], we 
consider the commuting set {J 2 ,v ra }, where the operator v ra is defined by 

v ra :=e^ r \j,-j)(j,j\+ Q U - m)a \j,m+l)(j,m\ (113) 

m=-j 



modulo its action on the space of constant angular momentum j spanned by the B^j+x 
basis. In f)113p . q is a primitive (2j + l)-th root of unity, i.e., 



q := e 2 ^( 2j+1 ) (114) 

and the parameters r and a are fixed parameters such that 

r eR, ae Z/(2j + 1)Z (115) 

It is to be noted that v ra is pseudo-invariant under the cyclic group CV,+i in the sense that 
it transforms as an IRC of C-y+i (different from the identity IRC). The common eigenstates 
of J 2 and v ra , associated with the SO (3) D C^j+i chain, provide an alternative basis to 
that given by the common eigenstates of J 2 and J z , associated with the SO (3) D SO (2) 
chain. This can be precised by the following result. 

Theorem 3. For fixed j, r and a, the 2j + 1 common eigenvectors of v ra and J 2 can 
be taken in the form 

i j 

\ja;ra) = = g (j+m)(j - m+1)a/2 ^ mr+(j+m)a |j, m) (116) 

* m=—j 

with a — 0, 1, . . . , 2j. The corresponding eigenvalues of v ra are given by 

v ra \ja;ra) = q^~ a \ja;ra) (117) 

so that the spectrum of v ra is non degenerate. 
The inner product 

(ja;ra\j/3;ra) = 6 a> p (118) 

shows that 

B ra :={\ja;ra):a = 0,l,...,2j} (119) 

is an orthonormal set which provides a nonstandard basis for the irreducible representation 
matrix of SU(2) associated with j. For fixed j, there exists a (2j + l)-multiple infinity of 
orthonormal bases B ra since r can have any real value and a, which belongs to the ring 
Z/ (2j + 1)Z, can take 2j + 1 distinct values (a = 0, 1, ... , 2j). 

4-3 Other bases in quantum information 

We now go back to quantum information. By using the change of notations 

d :— 2j + 1, n := j + m, \n) :— \j, — m), \aa; r) := \ja; ra) (120) 



adapted to quantum information and in agreement with (I109p . the operator v ra can be 
rewritten as 

d-l 



Vr 



e i ff (d-i)r| d _ i)(o| + ^g na | n -l}(n| 



[121) 



n=l 



Each of the eigenvectors 



\aa; r 



1 <Z-1 

) = g (d-l) 2 r-/4_L gn(d-»)tt/2+»[a-(d-l)r/2] | rf _ j _ ^ 



(122) 



n=0 



(with a = 0, 1, . . . , d — 1) of linear combination of the qudits |0), |1), . . . , \d — 1). 

For fixed d, r and a, the orthonormal basis 



B ra '■= {\aa; r) : a = 0, 1, . . . , d — 1} 



(123) 



is an alternative to the computational basis. As already mentioned, there is (i-multiple 
infinity of orthonormal bases B ra . 

All this can be transcribed in terms of matrices. Let V ra be the d x d matrix of 
the operator v ra . The unitary matrix V ra , builded on the basis Bd with the ordering 



0, 1, . . . , d — 1 for the lines and columns, reads 

/ q a 
q 2a 



y e m(d-l)r q q 
d-l 



Vr, 



The eigenvectors of V ra are 



Q 



\ 



(d-l)a 

o / 



(124) 



(aa; r) = g^" 1 ^ * V- g n( rf -n)a/2-n 



— n(d— l)r/2+na 



V-l-n 



(125) 



(with a = 0, 1, . . . , d — 1), where the <frk with k — 0, 1, . . . , d — 1 are the column vectors 

(126) 







W 



, • • • , 



-1 •= 



W 







representing the qudits |0), . . . , \d— 1), respectively. They satisfy the eigenvalue equa- 
tion 



V ra (p(aa; r) = q (<*-V(r+a)/2- a ^ aa . r ) 



(127) 



with a — 0, 1, . . . , d— 1. The Ka matrix can be diagonalized by means of the H ra unitary 
matrix of elements 



ITT \ 1 (d-l-n)(n+l)a/2+(d-l) 2 r/4+(d-l-n)\a-(d-l)r/2] 



(128) 



with the lines and columns of H ra arranged from left to right and from top to bottom in 
the order n, a — 0, 1, . . . , d — 1. Indeed, we have 



(H ra )^ V ra H ra = q 



(d-l)(r+o)/2 



fq° ... \ 

q- 1 ... 



\ 



(129) 



<d-l) J 



in agreement with (I127p . As an illustration, let us consider the d = 2 and d = 3 cases. 

For d = 2, we have two families of bases: the B t q family and the B r \ family (a can 
take the values a = and a = 1). The matrix (see (I124p ) 



Vrn 



q a 



q = e 



has the eigenvectors (see (I125P ) 

<t>(aa;r) = -y=(q' 



a/2-r/A+a 



o + g r/ », a = 0,1 



which correspond to the basis B ra . For r = 0, the bases 

Ax, : 0(00; 0) = (fa + O ) , 0(01; 0) = (0x - O ) 



(130) 



;m) 



(132) 



B i: 0(lO;O) = - ? =(0 1 + #o) 



;il;O) = - ? =(0 1 -# o ) 



(133) 



are (up to a rearrangement) familiar bases for qubits. 

For d = 3, we have three families of bases, that is to say B r0 , B r \ and B r2 , since a can 
be 0, 1 and 2. In this case, the matrix 

q a \ 
V,.„ := 10 q 2a 



q = e 



i2n/3 



/ 



admits the eigenvectors 



(aa: r 



r I a+2a— 2r 



r ( 

q (q 



>o + q a+a ~ r h + 2 ) , a = 0,1, 2 



For r = 0, the bases 



(134) 



(135) 



£? 00 : 0(OO;O) = - ? =(0 2 + 1 + 0o) 

0(Ol;O) = -^=(0 2 + g0 1 + g 2 0o) 
0(O2;O) = -^=(0 2 + g 2 0i + g0o) 



(136) 



Boi : 


0(10; 0) = 




02 + q<j>i + #o) 




0(11; 0) = 




02 + g 2 0i + <f>o) 




0(12; 0) = 




02 + 0i + g 2 0o) 




0(20; 0) = 




02 + g 2 0i + g 2 0o) 




0(21; 0) = 




02 + 01 + <?0o) 




0(22; 0) = 




02 + <?01 + 0o) 



(137) 



(138) 



are useful for qutrits. 

4-4 Mutually unbiased bases 

Going back to the case where d is arbitrary, we now examine an important property of 
the couple (B ra , B d ) and its generalization to couples (B ra , B rb ) with b ^ a. For fixed d, 
r and a, HI 22 j) gives 

Vra, a G {0, 1, . . . , d - 1} : | (n|aa; r)| = -= (139) 

Va 

Equation (I139P shows that B ra and -B^ are two unbiased bases. (Let us recall that two 
distinct orthonormal bases B a = {\aa) : a = 0,1,..., d— 1} and Bj, = {\bj3) : j3 = 
0, 1, . . . , d — 1} of the Hilbert space C d are said to be unbiased if and only if the inner 
product (aa\b(3) has a modulus independent of a and (3.) 

Other examples of unbiased bases can be obtained for d = 2 and 3. We easily verify 
that the bases B r0 and B r \ for d = 2 given by (11311) are unbiased. Similarly, the bases B r0 , 
B r \ and B r 2 for d = 3 given by f II 3 5 f) are mutually unbiased. Therefore, by combining 
these particular results with the general result implied by (I139P we end up with 3 mutually 
unbiased bases (MUBs) for d = 2 and 4 MUBs for d = 3. This is in agreement with the 
theorem according to which the number N MUB of pairwise MUBs in C d is such that 
3 < N MUB < d + 1 and that the maximum number d + 1 is attained when d is a prime 
number p or an integer power p e (e > 2) of a prime number p [1211 H22[ 1123] . The results 
for d = 2 and 3 can be generalized in the case where d is a prime number. This can be 
precised by the following theorem [TJl CC22 [E2SJ E23 fT2"B] . 

Theorem 4. For d = p, with p a prime number, the bases B r0 , B r i, . . . , B rp _i, B p 
corresponding to a fixed value of r form a complete set of p + 1 MUBs. The p 2 vectors 
\aa; r), with a, a; = 0, 1, ... ,p — 1, of the bases B r0 , B r i, . . . , -B rp _i are given by a single 
formula, namely (11221) or (11251) . The index r makes it possible to distinguish different 
complete sets of p + 1 MUBs. 



The proof is as follows. First, according to f)139p . the computational basis B v is 
unbiased with any of the p bases B r0 , B r i, . . . , B rp _i. Second, we get 

p-i 

(140) 



{aa: r 



1 P 

\bp;r) = ~^2q 



k(p-k)(b-a)/2+k(0-a) 



or 



V 



p-1 



k=0 



[aa;r\bP;r) = l^ e i7r{(a 

^ fc=0 



-b)k 2 + [p(b-a)+2(fS-a)]k}/p 



(141) 



The right-hand side of f)14ip can be expressed in terms of a generalized quadratic Gauss 
sum [129] 

(142) 



S(u,v,w) := e 



in(uk 2 +vk)/w 



k=0 



where u, v and w are integers such that u and w are mutually prime, uw is non vanishing 
and uw + v is even. This leads to 



with 



(aa;r|6/3;r) = -S(u,v,w) 
p 



u := a — b, v := —(a — b)p — 2(a — /3), w := p 



(143) 



(144) 



The generalized Gauss sum S(u,v,w) in (1143p - (ll44j) can be calculated from the methods 
described in |129j . We thus obtain 



\(aa;r\b/3;r}\ = — 



which completes the proof. □ 



(145) 



4-5 Mutually unbiased bases and Lie algebras 

4.5.1 Weyl pairs 

The matrix V ra can be decomposed as 



V ra — P r XZ a 



where 



P r :~- 



(\ 

1 

1 

\0 







iir(d— l)r 



(146) 



(147) 



and 



X :-- 



/O 1 
1 

OOO 
\1 





1 

0/ 



/l 

q 

q 2 

\o 



\ 




q d-l J 



(148) 



The linear operators corresponding to the matrices X and Z are known in quantum 
information as shift and clock operators, respectively. The unitary matrices X and Z 
g-commute in the sense that 



XZ - qZX = 



In addition, they satisfy 



X° 



(149) 



(150) 



where Id is the d- dimensional unit matrix. Equations ( I149P and (I150p show that X and 
Z constitute a Weyl pair |130j . The (X, Z) Weyl pair turns out to be an integrity basis 
for generating a set {X a Z b : a,b = 0, 1, . . . , d — 1} of d 2 generalized Pauli matrices in d 
dimensions (see for instance [1271 Q3H EM 033] in the context of MUBs and [TM fl35l fl36] 
in group-theoretical contexts). In addition, the set {q a X b Z c : a, b, c = 0, 1, . . . , d — 1} 
generates, with respect to matrix multiplication, a finite group of order d 3 , the Pauli 
group [127j . As an example, for d = 2 we have 



X = a x , Z = a z 



XZ 



-i(Ty, X°Z° = o"o 



(151) 



in terms of the ordinary Pauli matrices <j = I 2 , cr x , u y and cr 2 , and the Pauli group P 2 is 
isomorphic with the hyperbolic quaternion group. 

Equations (11491) and (I150p can be generalized through 



V ra Z - qZV ra = 0, e" 



-in(d-l)(r+a) 



(V ra ) d 



so that other pairs of Weyl can be obtained from V ra and Z. Note that 



X = Vqq, z={v 00 yvc 



01 



(152) 



(153) 



which shows a further interest of the matrix V rn - 



4.5.2 MUBS AND THE SPECIAL LINEAR GROUP 

In the case where d is a prime integer or a power of a prime integer, it is known that 
the set {X a Z b : a, b = 0, 1, . . . , d - 1} \ {X°Z } of cardinality d 2 - 1 can be partitioned 



into d + 1 subsets containing each d — 1 commuting matrices (cf. |131j ). Let us give an 
example. 

For d = 5, we have the 6 following sets of 4 commuting matrices 

V := {01,02,03,04}, Vi := {10, 20, 30, 40} 

V 2 := {11,22,33,44}, V 3 := {12, 24, 31, 43} (154) 

V 4 := {13,21,34,42}, V 5 := {14, 23, 32, 41} 

where ab is used as an abbreviation of X a Z b . 

More generally, for d = p with p prime, the p + 1 sets of p — 1 commuting matrices 
are easily seen to be 

V := {X°Z a :a=l,2,...,p-l} 
Vi := {X a Z°:a = l,2,...,p-l} 
V 2 := {X a Z a : a = 1,2,..., p- 1} 

V 3 := {X a Z 2a : a = 1,2,..., p- 1} (155) 

V p _! := {X a Z^- 2 ) a :a = l,2,...,p-l} 
V p := {X a Z ip - 1)a :a = l,2,...,p-l} 

Each of the p + 1 sets Vo, Vi, . . . , V p can be put in a one-to-one correspondence with one 
basis of the complete set of p + 1 MUBs. In fact, Vo is associated with the computational 
basis while V\, V 2 , ■ ■ ■ ,V P are associated with the p remaining MUBs in view of 

V 0a E V a+1 = {X b Z ab :b=l,2,...,p-l}, a = 0, 1, . . . ,p - 1 (156) 

Keeping into account the fact that the set {X a Z b : a, b = 0, 1, . . . ,p — 1} \ {X°Z } spans 
the Lie algebra of the special linear group SL(p, C), we have the following theorem. 

Theorem 5. For d = p, with p a prime integer, the Lie algebra sl(p, C) of the group 
SL{p, C) can be decomposed into a sum (vector space sum) of p + 1 abelian subalgebras 
each of dimension p — 1, i.e. 

sl(p, C) ^ v \£ v 1 \S . . . \£ v p (157) 

where the p + 1 subalgebras v ,Vi, . . . ,v p are Cartan subalgebras generated respectively 
by the sets Vq,V\, . . . ,V P containing each p — 1 commuting matrices. 

The latter result can be extended when d = p e with p a prime integer and e an integer 
(e > 2): there exists a decomposition of sl(p e ,C) into p e + 1 abelian subalgebras of 
dimension p e - 1 (cf. [Ml TM HSU ESS]). 



5 Appendix: The Racah parameters 



In the case of the I configuration, the Coulomb Hamiltonian Tic can be written as 

H c = (2£ + l) 2 £ Fk {[ n ^) 2 E( u(fc) «- u(fc) (j)) (158) 

k=0,2,...,2i ^ U U ^' i<j 

where the F k = Dk(£)Fk parameters are the usual Slater-Condon-Shortley parameters. 
It is clear that any linear transformation 

£ X = J2 b W^F k , X — 0,1, . . .,£ (159) 

fc=0,2,...,2l 

where b(£) is a regular matrix of dimension £+1 defines an equally acceptable parametriza- 
tion. 

As a trivial example, the D[. . .] parametrization in Section 2.5 corresponds to 
D[(00)0(M;)00] = (s\\u i0) \\s)-\£\\u^\\e)- 2 V2k + l(2£+lf ^ * o) F " (160) 

i.e., to a renormalization of the F k parameters. 

Less trivial examples are provided by the Racah parameters 

A = F° — -F 4 = F — 49F 4 
9 

B = 4il (9F2 ~ 5F4) = F2 ~ 5i?4 (161) 

C = — F 4 = 35F 4 
63 

for the d configuration |2] and the Racah parameters 

E° = F - 10F 2 - 33F 4 - 286F 6 

E l = i(70F 2 + 231F 4 + 2002F 6 ) 
9 

E 2 = l(F 2 -3F, + 7F 6 ) (162) 
y 

E 3 = \{5F 2 + 6F A - 91F 6 ) 

F — F° F — ! F 2 F — 1 t?4 ci 25 
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for the f N configuration [4J. The term energies for d N assume, to some extent, a simple 
form when expressed as functions of A, B and C . The E^ parameters (with j = 0, 1, 2, 3) 
for f N allow to decompose T-Lq m t° parts having well-defined properties under the action 
of the groups of the 50(7) D G 2 D 50(3) chain. 

As a last example, let us consider the parametrization defined by fl 1 5 9 [) with 

(I k £\ ( £ k £\ 
b(£) Xk = (-l)\2£+l)( n n )( x n ) (163) 



0/V-A A 



In this parametrization, the He operator can be rewritten as 

i 

u c =J2 v * ( 164 ) 

x=-i 

with 

V X = S X (2k + l)b(£) Xk J2(u {k) (i)-u {k) ti)) (165) 

k=Q,2,...,2l i<j 

We of course have V\ = V-\ and therefore there are I + 1 independent components V\ in 
Tic- The £ x parametrization was investigated in [1391 H40[ 1141] . Let us simply mention 
that the part Vq of He corresponds to a sum of surface delta interactions and that Tic 
can be reduced to Vq for 

F k = {2k + 1)F° (166) 

for k = 0, 2, . . . , 2£. In the special case of the d N configuration, it is to be realized that 
relation H166[) corresponds to the Laporte-Platt degeneracies |142] (see also [9| I139j H4"3] ) 
which occur for B = 0. 

6 Closing remarks 

Starting with the idea to substitute for the numerical methods of Slater, Condon and 
Shortley general methods close both to Dirac's ideas on quantum mechanics and to those 
of Wigner about the use of symmetries in physics, Racah developed practically in 20 years 
universal methods (irreducible tensor methods and group theoretical methods) used in 
many fields of physics and chemistry. In particular, the application of Racah's methods in 
atomic, nuclear and elementary particle physics as well as in group theory (Wigner-Racah 
algebra, state labeling problem) are well-known. We have shown how the use of Racah's 
methods in conjunction with SU(2) D G* or 50(3) D G symmetry adapted bases and 
effective operators yields sophisticated models in crystal- and ligand-field theories. In last 
analysis, these models are fully described by chains of groups, viz., the U(5) D SO (5) D 
SO (3) D G chain for the d N configuration in G and the U(7) D SO(7) D G 2 D SO(3) D G 
chain for the / configuration in G. 

As an application of current interest in the present days, we have shown the importance 
of the chain SO (3) D Cd for deriving a complete set of mutually unbiased bases when 
d is a prime integer. These bases are very useful in quantum information (quantum 
cryptography, quantum state tomography, quantum error codes) and equally in quantum 
mechanics (discrete Wigner function, mean King problem, path integral formalism). 

A common denominator to Sections 2, 3 and 4 is the notion of "chains of groups". 
Although chains of groups were in use before Racah (e.g., see the works by Ehlert on 



CH4 [144] . Bethe on crystal- field theory [29J and Wigner on supermultiplets of nuclei 
[145] ). his contribution to that part of applied group theory is essential and represents 
one of its major achievement so The interest for Physics and Chemistry of chains in- 
volving (noncompact and/or compact) continuous as well as finite groups is now well 
established. Such chains turn out to be useful in the investigation of broken symmetries 
which may arise either via descent in symmetry (Zeeman effect, homogeneous and inho- 
mogeneous Stark effect, ligand- field effect, etc.) or via spontaneous symmetry breaking 
(Landau and Jahn- Teller effect, symmetry breaking in elementary particle physics, etc.). 
In Racah's approach, which excludes the cases of external or Lorentzian and internal 
or gauge (super) symmetries, one group of the chain is a high symmetry group corre- 
sponding to a zeroth order approximation (like the cubic group in ligand-field theory) 
and another one is a low symmetry group corresponding to a first order approximation 
(like the tetragonal or trigonal group in ligand-field theory). The two symmetry groups 
correspond to known or postulated symmetries depending on whether the nature of the 
interactions involved is known or unknown. According to Wigner's theorem [147] . these 
symmetry or invariance groups (which leave invariant an Hamiltonian operator) provide 
representation labels or good quantum numbers for describing the state vectors. The 
other groups of the chain are dynamical or noninvariance groups in the sense that not 
all of its generators or elements commute with the Hamiltonian. They can describe part 
of the interactions and are generally introduced to make the chain as multiplicity-free 
as possible. Finally, the various groups of the chain are used to classify the state vec- 
tors and the (known or postulated) physical interactions. When elaborating a model 
based on symmetry considerations, the latter point is of considerable importance from 
a qualitative point of view (for level splitting and for selection rules) and a quantitative 
point of view (for the calculation of energy or mass matrices and transition probabilities). 
The preceding considerations apply to nuclear, atomic, molecular and condensed matter 
physics and also to quantum chemistry (chains of groups are even useful for classifying 
chemical elements |148j 1149] I15U] ). Note that the situation is a bit different in elementary 
particle physics since the notion of classification groups (with the pioneer works by Heisen- 
berg, Sakata, Gell-Mann Ne'eman and Zweig going from the SU(2) isospin group to the 
SU (3) D SU (2) chain involved in the first quark model) evolved to gauge groups (going 
from the SU (3)®SU {2)®U{1) D SU (3)®f7(l) standard model to the grand unified models 
based on the E s D E 7 D E 6 D 50(10) D SU(5) D SU(3)®SU(2)®U(1) D SU(3)®U(1)). 
However, in any field of physics there is a common scheme, namely, 



1 A fundamental result proved by Racah is that for a chain of groups having for head group a Lie group 
of order r and rank I, one can associate a complete set of commuting operators of cardinal (r + l)/2 (i.e., I 
Cartan operators plus / Casimir operators plus (r — 3Z)/2 labeling operators, some of the operators being 
Casimir, Cartan or labeling operators of the chain) [5]. See [146] for recent developments on this subject. 



(super) symmetries — > chain of groups — > invariance or 
co- variance — > conservation laws or good quantum numbers. 



To close this paper, let us add some further comments. Racah founded the main school 
of theoretical physics in Israel. He had a strong impact on (national and international) 
committees and on various research groups in theoretical and experimental spectroscopy 
(including the Laboratoire Aime Cotton in France). 

Racah had many students who deeply contributed to atomic and nuclear spectroscopy; 
they are profusely quoted in the review by Zeldes |151] . We would like to complete the list 
of students in the bibliography of Zeldes with a few words about Moshe Flato (Tel Aviv 
1937 - Paris 1998), a student of Racah during the period 1959-1963, who contributed to 
spread the ideas of Racah on crystal- and ligand-field theories. Flato achieved his M.Sc. 
thesis under the supervision of Racah in 1959 and prepared in 1960-1963 a Ph.D. thesis 
on a subject of nuclear physics (dealing with the Sp(2n) D U(n) chain in connection 
with the Elliott-Flowers model) given by RacahJ^I When Flato came to France in 1963 
his interest shifted to noncompact groups. He started working on Lorentzian symmetries 
and strong interactions when he was at Institute Henri Poincare in Paris (1963-1964). 
After that he was Associate Professor of Physics for three years at Universite de Lyon, 
then moved to Dijon in mathematics. Flato got in 1965 a Doctorat es Sciences Physiques 
from Universite de Paris on the basis of his works on elementary particle physicsjj Flato 
pursued a brilliant career both in France and worldwide, dealing with a great variety of 
subjects in physics and in mathematics. Flato evolved from theoretical physics to math- 
ematical physics and mathematics^ Among his many interests and contributions, let us 
mention the following: mass formulas (in relation with internal and external symmetries), 



2 According to Daniel Sternheimer, Flato was Racah's preferred student, probably the most brilliant 

in his generation. The families were friends since WWII when Flato's father was chief engineer of the 

British Mandate in Jerusalem. When Racah became Rector, he asked Flato to deliver (during 2 years, 

while Flato was doing military service) the traditional Racah lectures on group theory in physics, and 

recommended him for a course on solid state physics at Bar Ilan University, which Flato delivered without 

the compulsory yarmulke to students about his age. 

3 According to Sternheimer, after the death of Racah in 1965 (in Firenze on his way to join Flato in 

Paris) Flato decided not to publish the joint paper which they were preparing in nuclear physics. Racah 

had taken the manuscript with him in Firenze and intended to finalize it in Paris. Flato did not either 

defend in Jerusalem his Ph.D. based on that paper. Anyhow that became moot since he had already a 

French D.Sc. 

4 According to Sternheimer, his coworker for 35 years who heard with him a course on the theory of 
distributions by S. Agmon in Jerusalem in 1958-59, Flato had a dual training in physics and mathematics. 
Before opting for Racah he had considered working with S. Amitsur in algebra or with N. Rosen in 
relativity. 



conformal field theories; infinite-dimensional representations of Lie groups, singletons, 
AdS4/CFT 3 , composite electrodynamics; nonlinear representations of groups, covariant 
PDEs, global existence theorems for field theories (Yang-Mills, Maxwell-Dirac); and espe- 
cially the role of deformations in physics, including the now 35 years old and still frontier 
area of deformation quantization (symplectic and Poisson manifolds), quantum groups 
and noncommutative geometry. In 1968 Flato founded what was later called the Labora- 
toire Gevrey de Mathematique Physique at Universite de Bourgogne and in 1975 Letters 
in Mathematical Physics and two series of books published by Reidel (Kluwer). He had 
numerous students in France and abroad and a strong impact on Society (IAMP, Marie 
Curie chairs, Scientific Council of UAP). For more details, see |152[ 1153] . 

Racah and Flato shared important scientific and human qualities. Both were excellent 
teachers and at the same time exceptional researchers with a good sense of the duality 
theory-experiment, convinced of the importance of symmetries in physics. They knew how 
to communicate enthusiasm, give the right impulse to their students and collaborators, 
and inspire them to solve problems. Both had a strong impact on scientific communities 
and on national and international committees and enterprises. We learned and can still 
learn many things from them both from the human and scientific points of view. Their 
impact will last for a long time. We shall not forget them. 
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